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Abstract 

This  report  annotates  the  cold  regions  mobility  prediction  routines  included  in 
the  CAMMS/ALBE  mobility  models.  It  further  explains  the  development  of  the 
algorithms  that  are  used  in  these  models  to  describe  the  interaction  of  a  vehicle 
with  terrain  that  has  been  affected  by  cold  weather.  The  following  terrain 
canditions  are  discussed:  undisturbed  snow  (shallow  and  deep);  disturbed 
snow  (moderately  trafficked  and  hard  packed);  ice;  and  thawing  soils.  Several 
combinations  of  substrates  are  also  considered.  A  stand-alone  computer 
model  is  included. 


For  conversion  of  SI  metric  units  to  U.S./British  customary  units  ot  measurement 
consult  ASTM  Standard  E380-89a,  Standard  PracticeforUseof  the  International 
System  at  Units,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 
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Cold  Regions  Mobility  Models 

PAUL  W.  RICHMOND,  SALLY  A.  SHOOP  AND  GEORGE  L.  BLAISDELL 


INTRODUCTION 

Over  the  past  several  years  several  empirical 
and  theoretical  studies  have  aimed  to  predict  the 
mobility  of  vehicles  over  winter  terrain  and  the 
trafficability  of  that  terrain.  Most  of  these  efforts 
were  primarily  geared  toward  supporting  the  U.S. 
Army  Corps  of  Engineer's  AirLand  Battlefield  En¬ 
vironment  (ALBE)  demonstration  program. 
Many  of  these  models,  or  portions  thereof,  have 
also  been  proposed  for  incorporation  into  the 
NATO  Reference  Mobility  Model-II  (NRMM-II) 
( Ahlvin  and  Haley  1992)  and  the  Condensed  Army 
Mobility  Model  System  (CAMMS)  (Falls  et  al. 
1989).  These  general  mobility  models  can  be  used 
to  estimate  the  absolute  and  relative  mobility  of  ve¬ 
hicles  traveling  over  a  variety  of  terrain  types.  In 
addition  to  being  used  as  military  planning  and 
operational  tools,  the  models  can  also  be  used  to 
compare  specific  vehicles  or  to  determine  specific 
vehicle-terrain  interactions. 

In  this  report  the  interaction  between  winter 
terrain  surface  conditions  and  wheels  (driven  or 
undriven)  or  tracks  is  discussed  for  snow,  ice, 
freezing/ thawing  ground  and  layered  combina¬ 
tions  of  these  conditions.  This  report  brings  togeth¬ 
er  supporting  documentation,  data  and  theory 
for  individual  segments  of  the  cold  regions  mod¬ 
els.  A  stand-alone  program  and  the  FORTRAN 
code  for  CAMMS/NRMM  implementation  of  the 
cold  regions  models  are  included  as  appendices. 

BACKGROUND 

The  term  "mobihty"  is  defined  as  the  efficiency 
with  which  a  vehicle  travels  from  one  point  to  an¬ 
other.  Trafficability  refers  to  the  ability  of  the  ter¬ 
rain  to  support  mobility.  An  effective  mobility 
model  must  be  able  to  predict 

•  If  a  vehicle  can  propel  itself  between  two 
points  (go/no-go); 


•  If  the  vehicle  can  maneuver  or  tow'a  trailer 
(tractive  reserve);  and 

•  How  fast  the  vehicle  can  get  where  it's  go¬ 
ing  (speed  made  good). 

The  first  two  issues  can  be  determined  at  the 
same  time,  since  a  positive  value  for  tractive  re¬ 
serve  is  required  for  a  "go"  condition. 

In  its  simplest  form,  mobility  can  be  described 
by  a  single  equation  expressing  the  balance  be¬ 
tween  traction  T  and  resistance  R.  Tractive  re¬ 
serve  or  net  traction  can  be  calculated  from 

^net  ~  "^gross  ~  ^terrain  “  -^internal  ■  (1) 

In  eq  1,  is  the  maximum  tractive  force  that  a 
specific  vehicle  is  able  to  generate  on  a  particular 
terrain.  Gross  traction  is  generally  a  function  of 

•  The  tire/ track  contact  pressure; 

•  The  ability  of  the  running  gear  to  engage 
with  the  terrain  (e.g.  the  effect  of  tire  tread 
and  traction  aids,  or  track  grousers  and 
cleats); 

•  The  shear  strength  of  the  top  layers  of  the 
terrain;  and 

•  The  power  available  to  the  tire  or  track. 

Motion  resistance  can  be  divided  into  two  parts: 

that  produced  by  external  forces  and  that  pro¬ 
duced  by  internal  forces.  The  external  resistance 
^terrain  ^^e  resistance  attributable  to  the  surface 
and  is  a  function  of  the  strength  of  the  terrain  and 
the  vehicle's  running  gear  characteristics.  To  ob¬ 
tain  Rterrain  usually  necessary  to  determine 
the  level  to  which  a  vehicle  sinks  below  the  ter¬ 
rain  surface  (sinkage,  z).  Rintemal  is  the  resistance 
caused  by  friction  within  the  vehicle  (tire  defor¬ 
mation  or  track  roller  resistance,  friction  in  driv¬ 
eline  components,  etc.). 

Speed  made  good,  besides  depending  on  ter¬ 
rain  and  vehicle  characteristics  (slip  and  transmis¬ 
sion  power  curves),  is  also  strongly  dependent  on 
operator  skill,  visibility,  terrain  roughness  and 
other  parameters.  Although  important,  speed  made 
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Figure  1.  Determining  whether  a  snowpack  is  considered  shallow  or  deep. 


good  has  not  been  fully  investigated  under  win¬ 
ter  conditions.  Except  for  an  approximation  of 
speed  reduction  caused  by  vehicle  slip  in  snow 
(discussed  below),  it  is  not  addressed  by  any  spe¬ 
cific  cold  regions  models.  Effects  on  speed  deter¬ 
mined  for  temperate  climatic  conditions  should 
apply  to  winter  as  a  first  approximation. 

Vehicle  mobility  in  snow  depends  on  snow 
depth,  snow  density,  crystal  structure  of  the  snow 
grains,  liquid  water  content,  stratification  in  the 
snowpack  and  virtually  anything  that  affects  the 
mechanical  characteristics  of  the  snow  cover.  An 
overview  of  snow  characterization  measurements 
used  for  vehicle  mobility  and  snow  pavements  is 
given  in  Shoop  (1993b),  Shoop  and  Alger  (1993) 
and  Abele  (1990).  Snow  strength  characterization 
is  reviewed  in  Shapiro  et  al.  (in  prep.).  To  fully 
describe  a  snowpack  in  terms  of  its  physical  and 
mechanical  properties  is  time  consuming  and  re¬ 
quires  sophisticated  equipment  and  specialized 
techniques.  Since  the  models  described  in  this  re¬ 
port  were  designed  for  tactical  use,  it  was  unreal¬ 
istic  for  the  models  to  require  as  input  many  of 
the  traditional  scientific  measurements.  The  algo¬ 
rithms  we  describe  here  rely  solely  on  the  most 
basic  snow  parameters:  depth  and  density. 

SHALLOW  SNOW 

Eor  mobility  purposes,  snow  depth  can  be  cat¬ 
egorized  as  either  shallow  or  deep,  the  difference 
being  based  on  the  response  of  the  snowpack  to 
the  load  applied  by  the  vehicle.  If  the  region  of 
disturbed  snow  under  the  vehicle's  running  gear, 
called  the  pressure  bulb,  extends  to  the  ground  or 
pavement  under  the  snow  cover,  then  the  snow  is 
considered  shallow.  If  the  pressure  bulb  is  suspend¬ 
ed  in  the  snowpack,  a  deep  snow  condition  exists. 
An  additional  constraint  for  the  shallow  snow  case 
for  our  models  is  that  the  sinkage  of  the  vehicle 


be  less  than  the  vehicle  ground  clearance.  Inter¬ 
estingly,  whether  a  given  snowpack  thickness  is 
deep  or  shallow  depends  on  the  strength  of  the 
snow  and  the  characteristics  of  the  vehicle  to  be 
operated. 

Since  an  assessment  of  whether  a  given  snow 
condition  is  shallow  or  deep  requires  detailed 
knowledge  of  the  applied  load  and  snowpack  re¬ 
sponse,  we  suggest  two  guidelines  for  making 
this  determination.  These  guidelines  are  based  on 
the  amount  of  sinkage  calculated  in  relation  to  ei¬ 
ther  the  tire  radius  or  the  undercarriage  clearance 
as  illustrated  in  Figure  1. 

Undisturbed  snow  on 
a  firm  substrate 

We  describe  undisturbed  snow  as  snow  that 
has  been  driven  over  by  less  than  four  vehicles. 
We  chose  this  value  based  on  a  limited  number  of 
measurements  of  pressure  bulb  density,  where 
we  noted  that  it  takes  several  wheel  passes  over 
the  same  track  to  approach  the  critical  density  of 
snow.  (Critical  density  is  generally  accepted  to  be 
0.55  Mg/m^,  which  represents  the  maximum  den¬ 
sity  that  can  be  attained  by  snow  grain  rearrange¬ 
ment.  An  increase  in  density  beyond  this  requires 
deformation  of  individual  ice  grains.)  Further,  to 
be  called  undisturbed  snow  it  must  have  a  pre- 
vehicle-passage  density  of  less  than  0.55  Mg/ m^. 
This  is  almost  always  the  case  for  snow,  even  in 
polar  regions. 

If  the  terrain  underlying  the  snow  has  a  rating 
cone  index  (RCI)  (Shoop  1993a)  value  high  enough 
to  fully  support  the  vehicle  of  interest  with  no  sink- 
age,  mobility  is  assumed  to  be  unaffected  by  the 
substrate.  We  consider  an  RCI  value  of  100  or 
greater  to  be  adequate  to  resist  sinkage.  Nearly  all 
paved  surfaces  and  frozen  soils  have  an  RCI  greater 
than  100.  For  undisturbed  snow  on  a  firm  substrate, 
sinkage  will  only  occur  in  the  snow,  and  calcula¬ 
tion  of  traction  and  resistance  is  straightforward. 
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Figure  2.  Vehicle  traction  on  snow. 

For  shallow,  undisturbed  snow  lying  on  top  of 
a  strong  substrate,  we  calculate  gross  traction 
(kN)  for  one  driven  wheel  or  track  element  i 
(Blaisdell  et  al.  1990)  as 


Tgrossi  =  0-851 


(2) 


where  Nj  is  the  normal  stress  (kPa)  on  driving  ele¬ 
ment  i  and  Aj  is  that  element's  contact  area  (m^). 
The  total  gross  traction  for  a  vehicle  is  the  sum  of 
the  traction  obtained  for  each  driven  wheel  or 
track.  This  equation  is  based  on  data  generated 
from  a  wide  range  of  vehicle  types  and  sizes 
(both  wheeled  and  tracked,  with  contact  pres¬ 
sures  ranging  from  10  to  250  kPa)  (Fig.  2).  The 
correlation  coefficient  for  this  equation  (r^)  is 
0.97,  with  data  from  the  standard  military  ve¬ 
hicles  falling  within  -9  to  16%  of  the  predicted 
value. 

The  equation  for  motion  resistance  (N)  in  un¬ 
disturbed  snow  is 


^terrain  =  68.083  (pg  b 


(3) 


a.  Wheel.  b.  Track. 

Figure  3.  Snow  and  vehicle  characteristics. 
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where  pg  =  density  of  the  snow  immedi¬ 
ately  in  front  of  the  vehicle's 
miming  gear  (Mg/m^) 
a  =  length  of  the  tire  or  track  in 
contact  with  the  snow  (m) 
b  -  maximum  width  of  the  tire  or 
track  (m)  (Fig.  3). 

This  resistance  value  is  applied  to  up  to 
four  axles  or  one  track  set  on  the  vehicles  of 
interest  (including  any  towed  units).  For 
unknown  reasons,  predictions  of  resistance 
for  the  Small  Unit  Support  Vehicle  (SUSV) 
are  most  reliable  when  only  the  leading 
track  set  is  considered  (Richmond  et  al. 
1990).  Recent  analysis  by  Richmond  (in 
prep.)  indicates  that  only  the  first  three  ax¬ 
les  passing  through  undisturbed  snow  are  impor¬ 
tant;  we  apply  this  equation  to  the  first  four  axles 
since  the  HEMTT  with  four  axles  was  included 
in  development  of  eq  3.  The  data  forming  the  ba¬ 
sis  of  eq  3  are  displayed  in  Figure  4.  Although 
the  trend  is  strong,  the  value  is  low  (0.39),  and 
research  refining  our  understanding  of  vehicle 
motion  resistance  in  snow  is  continuing.* 

The  length  a  of  a  tire  or  track  contact  with 
snow  in  eq  3  requires  the  calculation  of  vehicle 
sinkage  z  (Fig.  3).  This  is  accomplished  by  the  fol¬ 
lowing  equation  for  sinkage  in  undisturbed  shal¬ 
low  snow; 


=  h 


pf ; 


(4) 


where  z, 


'max  =  maximum  predicted  sinkage  (cm) 
in  snow  occurring  under  the  tire  or 
track  with  maximum  contact  pres¬ 
sure  pmax 

h  =  initial  snow  depth  (cm) 

Po  and  pf  =  initial  and  theoretical  final  densities 
(Mg/m3),  respectively  (Fig.  5). 

The  final  density  pf  is  a  fimction  of  the  ap¬ 
plied  load.  In  temperate  regions  with  sea¬ 
sonal  snow,  the  following  pf  values  are  used 
as  estimates;  they  are  based  on  analysis  of 
field  sinkage  measurement,  as  illustrated  in 
Figure  6: 


*A  new  resistance  algorithm  presented  by  Rich¬ 
mond  (in  prep.)  for  wheeled  vehicles  in  snow  will  be 
incorporated  in  an  NRMM  version  newer  than  2.5.0. 
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Figure  4.  Vehicle  motion  resistance  in  shalloxv  snow. 


Figure  5.  Prediction  of  vehicle  sinkage  in  shallow 
snow. 


Figure  6.  Vehicle  sinkage  analysis.  (From  Richmond  et  al.  1990.) 
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The  sinkage  calculation  is  also  used  to  esti¬ 
mate  whether  the  snow  is  deep  or  shallow. 

Once  sinkage  z  is  determined,  the  value  of  a 
(contact  length)  for  wheeled  vehicles  is 

a  =  r  arccos  [(r-z)  /  r]  (5) 

where  r  is  the  tire  radius  (cm  or  in.).  For  tracked 
vehicles  the  equation  is 


analysis  that  is  beyond  the  scope  of  our  initial 
models.  Further  research  may  show  that  these 
shortcomings  can  be  accounted  for  by  lowering 
the  RCI  value  used  to  separate  firm  from  soft  sub¬ 
strates.  We  have  also  ignored  the  possibility  that 
shear  zones  could  form  between  the  snow  and 
the  soil  (most  likely  in  shallow  snow),  which 
would  affect  both  traction  and  resistance  calcula¬ 
tions. 

Soil  sinkage  is  calculated  separately  for  wheeled 
and  tracked  vehicles  using  the  sinkage  equations 
of  Willoughby.  The  equations  for  wheeled  and 
tracked  vehicles,  respectively,  are 


a  =  z/sin  (0)  (6) 

where  6  is  the  angle  between  the  leading  edge  of 
the  track  and  the  ground  surface  (Fig.  3b).  Since 
the  vehicle  database  defined  for  CAMMS  and 
NRMM  does  not  contain  this  value,  we  assume 
an  average  value  of  26°  for  0  for  all  tracked  ve¬ 
hicles. 

Undisturbed  snow  on 
a  soft  substrate 

For  undisturbed  snow  over  a  soft  soil,  addi¬ 
tional  sinkage  occurs  as  a  soil  rut  is  formed.  This 
increases  motion  resistance  as  the  wheel  or  track 
sinks  deeper  in  the  snow  and  possibly  below  the 
soil  surface. 

In  this  case  the  total  sinkage  is  estimated  to  be 
the  sum  of  the  sinkage  of  the  two  materials  inde¬ 
pendently.  Soil  sinkage  is  calculated  ignoring  the 
snow  cover  and  is  added  to  the  sinkage  deter¬ 
mined  for  the  snow  from  eq  4.  The  arc  length  a, 
and  hence  the  motion  resistance  Rferraiiv  is  calcu¬ 
lated  based  on  the  combined  sinkage: 

^  “  ^soil  ^max-  (7) 

This  approach  may  be  naive  in  its  assumption 
that  the  snow  and  soil  can  each  be  treated  as 
separate  entities  and  that  their  effects  can  be 
added  to  determine  the  combined  effect  on  a  ve¬ 
hicle.  However,  this  assumption  is  probably  ad¬ 
equate  for  snow  depths  greater  than  about  15  cm 
and  soil  strengths  great  enough  to  suffer  less  than 
about  10  cm  of  sinkage.  We  also  recognize  that  we 
have  not  accounted  for  the  ability  of  snow  to 
spread  and  attenuate  the  vehicle's  load.  In  calcu¬ 
lating  soil  sinkage  using  this  approach,  we  have 
placed  the  entire  running  gear  load  and  contact 
pressure  on  the  soil.  Although  this  is  clearly  in¬ 
correct,  an  accurate  determination  of  the  actual 
load  and  pressure  would  involve  a  sophisticated 
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Zsoii  =  0.00443  Li  e(5  889  Wj /RCI  bj  Lp  (Tracked)  (9) 
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Zsoil  =  0  (If  RCI  >  100)  (10) 


where  RCI 
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rating  cone  index 

track  length  (in.) 

undeflected  tire  width  or  track 

width  (in.) 

tire  radius  (in.) 

wheel  or  track  load  (lb) 

tire  loaded  deflection  (in.) 

tire  unloaded  section  height  (in.) 

slip  in  decimal  form  (chosen  to  be 

0.05  for  Cold  Regions  Models)  for 

the  running  gear  element  i. 


Combined  sinkage  less  than  or 
equal  to  the  snow  depth 

For  the  situation  where  sinkage  occurs  in  the 
soil  (RCI  <  100)  but  z  is  still  less  than  the  snow 
depth  h,  we  assume  that  the  motion  resistance  is 
only  that  due  to  snow.  While  the  combined  z  from 
eq  7  is  larger  than  from  the  snow,  the  tire  or 
track  is  still  contained  completely  in  the  snow- 
pack  and  thus  its  resistance  to  motion  is  caused 
primarily  by  snow  deformation.  Equations  2  and 
3  are  used  to  calculate  traction  and  resistance. 


Personal  communication  with  W.  Willoughby,  U.S. 
Army  Waterways  Experiment  Station,  1992. 
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Figure  7.  Vehicle  traction  for  snow  over  ice.  (From  Richmond  et  al. 
1990.) 


Combined  sinkage  greater  than 
the  snow  depth 

For  some  situations  it  is  possible  that  the  com¬ 
bined  sinkage  z  is  greater  than  the  snow  depth.  In 
this  case  the  motion  resistance  due  to  the  sinkage 
into  the  soil  must  be  added  to  that  from  the  snow 
to  determine  the  total  resistance.  Unfortunately 
calculating  the  resistance  for  a  sinkage  greater 
than  the  snow  depth  involves  knowing  the  inter¬ 
action  of  the  snow  and  the  substrate  during  de¬ 
formation  and  shearing  and  is  far  beyond  the 
scope  of  current  understanding.  Thus,  we  have 
no  provisions  for  calculating  the  additional  resis¬ 
tance  or  the  effect  on  traction  when  z  is  greater 
than  h.  The  model  calculates  the  resistance  based 
on  the  total  sinkage  considered  to  be  in  snow,  and 
traction  is  based  on  the  condition  specified  as  if 
this  case  did  not  occur. 

Undisturbed  snow  over  ice 

The  algorithms  above  have  all  assumed  that 
shear  displacement  in  the  terrain  as  a  result  of  the 
tractive  demand  of  the  running  gear  occurs  with¬ 
in  the  snowpack.  In  this  assumption  we  are 
claiming  that  the  internal  shear  strength  of  the 
snow  below  a  tractive  element  is  less  than  or 
equal  to  the  shear  strength  of  the  substrate  or  the 
interface  between  the  substrate  and  the  snow. 
This  is  most  likely  true  for  all  situations  except 
when  the  snow  overlies  ice.  This  could  occur  on  a 
frozen  river  or  lake  or  on  a  road  surface  where 
precipitation  started  as  freezing  rain  prior  to  be¬ 
coming  snow.  In  this  case  the  interfacial  shear 
strength  between  the  snow  and  ice  is  almost 
guaranteed  to  be  less  than  the  shear  strength  within 


the  snow.  Thus,  for  shallow  snow  where  the  pres¬ 
sure  bulb  is  in  contact  with  the  ice,  we  assume 
that  traction  will  be  governed  by  the  interfacial 
properties. 

For  undisturbed  shallow  snow  deeper  than  1 
cm  lying  on  top  of  ice,  we  use  the  following  trac¬ 
tion  equation: 

Tg,oss,  =0.127^106^  (11) 

which  was  presented  by  Richmond  et  al.  (1990) 
(Fig.  7).  The  value  of  1  cm  was  chosen  to  reflect  the 
point  at  which  the  snow/ice  interface  would  have 
effectively  no  strength,  and  the  traction  generated 
would  be  the  same  as  that  on  a  clear  ice  surface, 
^terrain  ^his  case  is  calculated  using  eq  3. 

Disturbed,  processed 
and  packed  snow 

The  physical  and  mechanical  properties  of  a 
snowpack  change  significantly  when  it  is  driven 
over  repeatedly.  After  a  large  volume  of  traffic, 
the  snow's  characteristics  cease  to  change  dra¬ 
matically  as  the  result  of  tire  loads,  and  the  snow 
is  then  considered  to  be  packed. 

Vehicle  operators  tend  to  follow  in  the  tracks 
of  preceding  vehicles.  We  assume  that  after  four 
passes  in  the  same  track,  the  snow  has  reached  its 
critical  density,  and  succeeding  vehicles  will  be 
traveling  on  packed  snow.  Usually,  packed  and 
disturbed  snow  conditions  can  only  be  formed  over 
a  firm  substrate. 

In  some  cases  an  area  of  snow  is  mechanically 
processed  to  produce  a  snow  road  that  is  capable 
of  supporting  wheeled  vehicle  traffic  (Abele  1990). 
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Figure  8.  Vehicle  traction  for  packed  snow.  (From  Richmond  et  al. 
1990.) 


We  assume  that  disturbed,  packed  or  processed 
srrow  attains  a  density  of  at  least  0.55  Mg/ m^. 
Since  only  vehicles  that  have  a  very  high  ground 
pressure  are  expected  to  cause  any  sinkage  on 
these  surfaces,  Rterrain  this  case  becomes  zero 
for  all  vehicles. 

The  equation  we  use  to  predict  traction  on 
packed  snow  is 

Tgrossi  =  0.321N°-^^Ai  (12) 

and  is  shown  in  Figure  8. 

A  comparison  of  measured  and  predicted  trac¬ 
tion  on  shallow,  groomed  snow  is  discussed  in 
Appendix  C. 

For  disturbed  snow  deeper  than  1  cm  over  an 
ice  surface,  we  assume  that  a  bond  has  developed 
between  the  packed  snow  and  ice;  we  use  eq  12 
for  traction  in  this  case  and  again  assume  that  the 
Rsnow  zero. 

DEEP  SNOW 

Deep  snow  occurs  when  the  deformation  bulb 
under  a  wheel  or  track  does  not  extend  down 
through  the  snowpack  to  an  underlying  surface 
(Fig.  1).  Since  vehicle  sinkage  depends  on  ground 
pressure,  it  is  apparent  that  the  presence  of  a  deep 
snow  condition  is  a  function  of  the  snow  depth, 
bearing  capacity  (generally  reported  in  the  form 
of  a  pressure-sinkage  relationship)  and  vehicle 
characteristics.  Additionally  we  include  as  deep 
snow  the  cases  where  sinkage  is  greater  than  2/3 


of  the  wheel  radius  and  where  sinkage  is  greater 
than  the  ground  clearance  for  wheeled  or  tracked 
vehicles.  This  added  set  of  criteria  for  defining  a 
deep  snow  condition  recognizes  that  when  a  ve¬ 
hicle  sinks  to  the  point  where  the  majority  of  its 
running  gear  is  below  the  snow  surface,  added 
sources  of  resistance  arise  beyond  simple  com¬ 
paction  of  the  snow  in  front  of  the  tire  or  track. 

Both  traction  and  resistance  are  different  for 
deep  snow  than  for  shallow  snow.  However, 
there  are  few  published  results  for  deep  snow 
mobility  with  wheeled  and  tracked  vehicles, 
making  it  difficult  to  write  accurate  prediction 
equations.  Thus,  we  used  the  shallow  snow  re¬ 
sults  as  a  starting  point  for  the  deep  snow  mod¬ 
els.  We  began  by  assuming  that  the  deep  snow 
condition  only  alters  traction  to  a  small  degree, 
since  the  majority  of  the  tractive  force  is  gener¬ 
ated  in  the  fully  compacted  portion  of  the  pres¬ 
sure  bulb.  The  tractive  force  developed  in  deep 
snow  is  most  likely  somewhat  less  than  for  shal¬ 
low  snow,  because  the  pressure  bulb  will  have  a 
slightly  lower  density.  The  degree  of  reduction  in 
density  and  tractive  force  is  unknown,  so  we 
chose  to  use  the  shallow  snow  traction  expression 
(eq  2)  for  deep  snow. 

Motion  resistance  generated  as  the  result  of 
sinkage  into  the  snow  exists  in  deep  snow  just  as 
it  did  for  shallow  snow.  In  the  vast  majority  of 
cases,  sinkage  is  greater  in  deep  snow,  so  resis¬ 
tance  is  higher,  since  a  greater  portion  of  the  wheel 
or  track  is  involved  in  compacting  and  displacing 
snow.  However,  several  additional  sources  of  resis¬ 
tance  may  arise  in  deep  snow.  For  example,  por- 
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For  wheeled  vehicles: 


Light  snow  over  smooth  ice  (Shoop  1993b) 
"Rough"  ice  (Shoop  1993b) 

“Rough"  ice  (Biaisdell  1993b) 

•  Buffed  ice  at  -3°C  (Biaisdell  &  Borland  1 992) 

•  Buffed  ice  at  - 1 0°C  (Biaisdell  &  Borland  1 992) 

•  Pi^^(SME  1993) 

•  "Rule  of  thumb" 

Highway  travel  (Clark  1982) 

I _ ^ ^ _ I  I 

0  0.1  0.2  0.3  0.4 

Traction  Coefficient  for  Ice 

Figure  9.  Vehicle  traction  on  ice. 

tions  of  the  undercarriage  may  drag  in  the  snow 
as  the  result  of  greater  sinkage.  In  very  deep  snow, 
even  the  front  face  of  a  vehicle  may  be  "plow¬ 
ing"  snow. 

Lacking  quantitative  data  for  motion  resis¬ 
tance  in  deep  snow,  we  based  our  initial  model  of 
resistance  in  deep  snow  on  a  combination  of  ana¬ 
lytical  considerations  and  our  collective  obser¬ 
vations  of  the  behavior  of  vehicles  operating  in 
deep  snow  conditions.  Our  model  of  resistance 
in  deep  snow  uses  "multipliers"  applied  to  the 
shallow  snow  resistance  expression  (eq  3).  These 
multipliers  were  chosen  to  represent  the  signifi¬ 
cant  increase  in  resistance  experienced  when  snow 
is  being  pushed  forward  by  the  near-vertical  face 
of  the  tire  when  sinkage  is  large.  For  tracked  ve¬ 
hicles  the  concern  is  primarily  interference  with  the 
undercarriage.  The  multipliers  are  as  follows: 


Po  <  0.15  Mg/m^ 

1)  If  sinkage  >  wheel  radius,  multiply 
resistance  by  1.5. 

2)  If  sinkage  >  ground  clearance,  add 
1.5  to  multiplier  above. 

po^  0.15  Mg/m^ 

1)  If  sinkage  >  (0.666  x  wheel  radius) 
but  less  than  wheel  radius,  multiply 
resistance  by  1.5. 

2)  If  sinkage  >  wheel  radius,  multiply 
resistance  by  2.5. 

3)  If  sinkage  >  ground  clearance,  add 
1.5  to  multipliers  above. 

For  tracked  vehicles: 

Po  <  0.15  Mg/m^ 

1)  If  sinkage  >  ground  clearance,  multi¬ 
ply  resistance  by  2.5. 

Po  ^  0.15  Mg/m^ 

1)  If  sinkage  >  ground  clearance,  multi¬ 
ply  resistance  by  4.0. 

In  formulating  these  guidelines,  we  took  ground 
clearance  to  be  the  height  above  ground  of  the 
lowest  major  component  of  the  undercarriage  of 
a  vehicle,  measured  while  it  is  parked  on  a  firm 
surface.  Our  experience  suggests  that  small,  or  a 
low  number  of,  protuberances  generate  little  ad¬ 
ditional  resistance  (Richmond,  in  prep.).  How¬ 
ever,  the  NRMM/ CAMMS  database  lists  ground 
clearance  as  the  distance  from  the  lowest  compo¬ 
nent  to  the  ground.  Therefore,  our  model  will 
yield  conservative  output  (it  will  err  on  the  side 
of  overprediction  for  Rferrain)- 


Figure  10.  Vehicle  traction  on  snow 
and  ice. 


Normal  Stress  (kPa) 
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ICE 

Ice  traction  is  a  function  of  surface  roughness, 
temperature,  tire  compound  and  vehicle  speed 
(Blaisdell  and  Borland  1992).  While  it  may  be  pos¬ 
sible  to  easily  incorporate  ice  temperature  into  the 
NRMM/CAMMS  database,  it  is  unlikely  that  ei¬ 
ther  tread  compound  or  surface  roughness  could 
be  included  with  any  degree  of  accuracy  in  the  near 
term.  So,  to  treat  the  case  of  traction  on  ice,  we  have 
reviewed  published  data  (Fig.  9)  (Shoop  1993b)  and 
taken  a  conservative  value  of  0.1.  Therefore,  we  cal¬ 
culate  traction  on  ice  as 

Tg,oss,  =  0.1NiAi.  (13) 

The  traction  equations  for  ice  and  snow  are  com¬ 
pared  in  Figure  10.  A  model  of  ice  traction  utiliz¬ 
ing  traction  aids  exists  (Blaisdell  1984),  but  it  re¬ 
quires  as  input  geometric  features  of  the  devices, 
which  are  not  available  in  the  current  Army  mo¬ 
bility  database. 

Because  of  the  essentially  undeformable  na¬ 
ture  of  ice  under  a  vehicle's  running  gear,  motion 
resistance  ^terrain  is  assumed  to  be  zero. 

FREEZING  OR  THAWING 
GROUND  CONDITIONS 

Freezing  ground  can  often  increase  vehicle  mo¬ 
bility,  while  thawing  ground  nearly  always  reduces 
mobility.  An  additional  issue  of  importance  is  the 
possibility  of  severe  terrain  damage  when  vehicles 
operate  in  areas  with  thawing  conditions.  Three 
critical  conditions  for  vehicle  mobility  on  freezing 
and  thawing  soils  are  illustrated  in  Figure  11. 


When  a  wet,  weak  soil  begins  to  freeze  (Fig. 

lla) ,  the  strong  frozen  layer  on  top  will  increase 
the  bearing  capacity  of  the  ground.  The  bearing 
capacity — the  ability  of  the  soil  to  support  a  ve¬ 
hicle — can  be  expressed  as  a  function  of  the  frost 
depth  t  and  the  relative  wetness  of  the  frozen 
ground  surface. 

During  spring  and  intermittent  thaws,  a  thawed 
layer  of  soil  develops  over  the  frozen  soil  (Fig. 

llb) .  In  the  top  layer  the  soil  moisture  is  higher 
than  normal  due  to  snowmelt,  rain  and  the  in¬ 
creased  water  drawn  to  the  soil  layer  during  the 
freezing  process.  This  moisture  is  trapped  in  the 
shallow  thawed  layer,  creating  a  wet  and  weak 
layer  of  soil  over  the  stronger  frozen  layer.  The 
reduction  in  vehicle  mobility  will  be  a  function  of 
the  strength  of  the  composite  soil,  which  can  be 
expressed  in  terms  of  thaw  depth  S  and  the  soil 
moisture  in  the  thawed  layer.  As  the  thaw  pro¬ 
gresses  deeper  (Fig.  11c),  the  frozen  layer  becomes 
too  distant  to  add  support  to  the  vehicle  or  strength 
to  the  effective  soil  system,  but  it  often  continues 
to  restrict  the  soil  drainage. 

Our  current  state  of  development  of  mobility 
modeling  for  freezing  ground  is  limited  to  go/ 
no-go  predictions  based  on  whether  the  ground 
can  support  the  vehicle.  For  thawing  ground  we 
can  make  more  quantitative  predictions  of  the 
effect  on  traction  and  motion  resistance.  Our 
models  assume  a  baseline  traction  and  motion 
resistance  for  the  soil  of  interest  in  its  fully 
thawed  state.  The  effects  of  freezing  and  thaw¬ 
ing  are  expressed  as  modifications  (multipliers) 
to  the  baseline  values  for  traction  and  motion 
resistance.  The  effects  of  other  factors  (e.g.  veg¬ 
etation  and  slope)  on  mobility  during  freezing 


a.  Critical  depth  of  frozen  soil  that  will  support  a 
vehicle. 


b.  Critical  depth  of  a  thawed  wet  layer  where  traction 
is  too  low.  If  the  tires  can  engage  the  strength  of  the 
frozen  layer,  then  the  ground  may  be  trafficable. 


c.  Frozen  layer  too  far  down  to  give  support.  The  layer 
still  impedes  drainage.  Moisture  content  and  soil  proper¬ 
ties  are  critical  for  trafficability. 


Frozen  Soil 


Unfrozen  Soil 


Figure  11.  Critical  conditions  for  trafficability  of 
freezing  or  thawing  ground. 
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or  thawing  are  assumed  to  be  the  same  as  for  un¬ 
frozen  conditions. 

Bearing  capacity  of 
freezing  ground 

Vehicle  operation  on  freezing  ground  is  char¬ 
acterized  by  the  ability  of  the  ground  to  fully  sup¬ 
port  the  vehicle.  For  unfrozen  soils  that  are  not 
extremely  difficult  to  traverse,  the  presence  of  5 
cm  of  frost  at  the  surface  will  usually  allow  un¬ 
limited  cross  country  operation  (Richmond  1991). 
Terrain  that  is  normally  untrafficable  may  require 
a  substantial  frost  layer  before  vehicle  operations 
are  possible. 

Certainly  "soft"  terrain  stands  to  gain  consid¬ 
erable  strength  upon  freezing.  For  example,  the 
compressive  strength  of  frozen  peat  can  be  350- 
400  times  its  unfrozen  strength  (MacFarlane  1968), 
making  it  possible  to  operate  vehicles  in  peat  areas 
as  long  as  they  are  frozen.  We  use  a  definition  of 
soft  terrain  based  on  the  Cone  Index,  where  val¬ 
ues  of  50  or  less  constitute  very  weak  soils.  Al¬ 
though  cone  penetration  is  not  particularly  appli¬ 
cable  to  measuring  the  strength  of  heavily  veg¬ 
etated  terrain,  it  can  serve  as  a  gross  assessment 
of  the  bearing  capacity  of  the  ground,  disregard¬ 
ing  the  strength  of  root  systems.  Niemi  and  Bayer 
(1970)  have  documented  Cone  Index  values  of  50 
or  less  for  peat  or  muskeg. 

To  arrive  at  predictions  of  bearing  capacity  for 
frozen  ground  over  a  soft  substrate,  we  utilize  pub¬ 
lished  guidelines  for  forestry  operations  on  peat- 
lands  (Rummukainen  1984)  (Table  1).  The  equiva¬ 
lent  loads  are  estimates  of  the  gross  vehicle  weight. 
From  these  guidelines  the  following  equations 
were  generated  to  represent  the  limits  for  a  break¬ 


through  failure  of  a  frozen  layer  that  is  less  than 
0.5  m  thick: 

P  =  0.35  for  dry  conditions  (14) 

P  =  0.86  for  wet  conditions  (15) 

where  P  is  the  maximum  recommended  load  in 
MN  and  t  is  the  frost  depth  in  meters.  The  wet 
condition  represents  saturated  peat.  In  the  dry  con- 


p 


a.  Tensile  failure  in  the  lower  beam  fibers. 


p 


f  1 


b.  Punching  through  the  frozen  layers. 


p 


c.  Crushing  and  compaction  of  the  frozen 
material. 

p 


d.  General  shear  failure  of  the  underlying 
material  and  tensile  failure  of  the  frozen 
layer. 


Figure  13.  Bearing  failure  modes  for  a 
shallow  frozen  layer  over  soft  ground. 
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Table  1.  Bearing  strength  of  frozen  peatland.  (After  Rummakainen  1984.) 


Thickness  of  frozen  peat  layer  (m) 
Dry  top  Wet  top 

peat  layer  peat  layer 

Bearing  capacity 

Approx,  equiv. 

load 

(MN) 

(Ibf) 

0.10 

0.05 

Will  bear  a  horse 

0.00445-0.008  1000-1800 

0.15-0.20 

0.10 

Will  bear  6-ton 

0.05 

12,000 

horse-drawn  sled  traffic 

0.20-0.35 

0.15-0.25 

Will  bear  empty 

0.036 

8000 

4-ton  truck 

0.35-0.50 

0.25-0.40 

Will  bear  10-ton  truck  traffic 

0.09 

20,000 

dition  the  peat  is  less  than  saturated  and  has  usu¬ 
ally  been  drained.  (A  frost  depth  of  0.5  m  will  pre¬ 
vent  breakthrough  failure  of  most  vehicles.)  Equa¬ 
tions  14  and  15  are  shown  schematically  in  Figure 
12  in  terms  of  the  frost  thickness  necessary  to  sup¬ 
port  different  vehicle  classes  (vehicle  class  is  ap¬ 
proximately  the  vehicle  weight  in  tons). 

Mobility  of  heavy  equipment  may  be  limited 
by  localized  crushing  failure  of  the  surface  mate¬ 
rial,  particularly  for  the  dry  conditions  (Fig.  13) 
(Shoop,  in  prep.).  Therefore,  the  prediction  equa¬ 
tion  for  dry  conditions,  where  surface  crushing  is 
likely,  is  not  applicable  for  vehicles  with  high 
ground  pressure.  We  recommend  that  the  dry  equa¬ 
tion  be  used  for  all  tracked  vehicles  and  wheeled 
vehicles  weighing  less  than  12  tons. 

The  predictive  formulas  are  based  on  a  best-fit 
equation  for  the  wheeled  vehicles,  and  they  over¬ 
estimate  the  frost  thickness  necessary  for  the  sled 
(Table  1).  Thus,  our  load  support  predictions  will 
be  conservative  for  vehicles  with  running  gear 
that  distributes  load  uniformly  (e.g.  tensioned 
tracks,  numerous  road  wheels,  skis). 

These  equations  can  also  be  applied  to  frozen 
ground  other  than  peat.  Once  frozen,  the  strength 
of  ground  depends  primarily  on  its  ice  content, 
its  density  and  its  temperature.  The  time-depen- 
dent  compressive  strength  of  frozen  peat  is  simi¬ 
lar  to  that  of  frozen  mineral  soils  (MacFarlane 
1968)  and  falls  within  the  range  of  frozen  silt  or 
clay  (Shoop,  in  prep.). 

Very  few  data  exist  for  vehicle  traction  on  frozen 
groimd.  In  general,  frozen  ground  enhances  mobil¬ 
ity  and  traction.  Exceptions  are  if  the  ground  has  a 
very  high  ice  content  or  if  the  temperature  is  near 
melting.  In  these  cases  traction  may  be  reduced  be¬ 
cause  of  the  slipperiness  of  the  surface  and  may  be 
closer  to  the  level  of  traction  experienced  on  ice 


rather  than  that  for  the  soil  alone.  However,  since 
we  have  no  data  to  confirm  this,  no  algorithm  is 
implemented  to  accoimt  for  this  effect. 

Effect  of  thawing  conditions  on 
vehicle  performance 

Thawing  ground  causes  problems  for  vehicle 
mobility  when  it  is  associated  with  thaw  weaken¬ 
ing  of  the  soil.  Ehiring  thaw,  otherwise  freely  drain¬ 
ing  soils  can  become  saturated  because  drainage 
is  reduced  by  the  underlying,  nearly  imperme¬ 
able  frozen  layer.  Vehicle  travel  is  then  restricted 
or  impossible,  and  traffic  can  cause  environmen¬ 
tal  damage  (torn  vegetation,  mass  soil  flow  and 
rutting,  and  subsequent  erosion  of  sediments). 

To  predict  traction  and  resistance  in  thawing 
soils,  we  start  from  either  measured  or  predicted 
vehicle  performance  for  the  same  soil  at  a  tem¬ 
perate  state.  Adjustments  are  made  to  compen¬ 
sate  for  the  significant  loss  of  shear  strength  in  a 
wet,  thawing  soil  (resulting  in  a  loss  of  traction) 
and  the  greater  vehicle  sinkage  in  the  thaw-weak¬ 
ened  soil  (causing  an  increase  in  motion  resistance). 
Using  established  methods  for  predicting  traction 
and  motion  resistance  for  unfrozen  soil  (in  this  case, 
models  developed  by  WES),  we  apply  multipliers 
to  reduce  traction  and  increase  motion  resistance 
as  the  result  of  the  thawing  condition: 

^thaw  ~  /  ^gross 

^thaw  “  S  ^terrain 

where  Tj^aw  =  gross  traction  in  the  thawing  soil 
Kthaw  =  terrain  motion  resistance  in  the 
thawing  soil 

/  =  traction  reduction  multiplier 
g  =  motion  resistance  multiplier 
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Traction  Coefficient 


1 


Figure  14.  Traction  coefficient  and  resistance  coefficient  vs.  water  content.  Below  the  liquid  limit,  traction  is  a  function  of 
applied  stress  as  reflected  in  tire  inflation  pressure.  The  net  tractive  coefficient  decreases  and  the  motion  resistance  coefficient 
increases  for  soil  water  contents  above  the  liquid  limit. 


Rterrairt  =  terrain  resistance  for  the  unfrozen 
soil. 

The  multipliers  are  defined  based  on  the  mea¬ 
surements  of  vehicle  performance  on  thawing  soil 
reported  in  Shoop  (1990, 1993c).  The  experiments 
were  performed  using  an  instrumented  vehicle  to 
measure  traction  and  motion  resistance  on  a  wide 
variety  of  thawing  conditions  for  a  fine-grained 
sand.  The  experimental  variables  were  soil  thaw 
depth,  moisture  content,  density  and  applied  stress 
(tire  inflation  pressure). 

The  multipliers  /  and  g  vary  with  soil  type, 
wetness,  contact  stress  and  thaw  depth.  Traction 
and  motion  resistance  vary  with  the  wetness  of 
the  thawed  soil  layer  as  shown  in  Figure  14.  At 
lower  water  contents  the  traction  coefficient  de¬ 
pends  on  the  tire  contact  stress.  At  high  soil  mois¬ 
ture  (above  the  liquid  limit),  traction  drops  rap¬ 
idly,  and  the  contact  stress  has  no  noticeable  ef¬ 
fect.  Motion  resistance  increases  considerably  at 
soil  moistures  above  the  liquid  limit,  and  the  ef¬ 
fect  of  tire  inflation  pressure  was  negligible  for 
the  conditions  tested. 

The  effect  of  the  depth  of  thaw  can  be  seen 
clearly  when  the  data  are  grouped  according  to 
the  wetness  of  the  thawed  layer  (Fig.  15).  The  de¬ 
gree  of  wetness  in  the  thawed  layer  is  defined  in 
Table  2.  For  both  traction  and  motion  resistance, 
the  "wet"  conditions  are  the  most  critical.  When 
the  soil  is  "moist,"  the  soil  strength  is  near  opti¬ 
mum,  and  much  of  its  strength  is  retained  during 
the  thaw  cycle.  Thus,  for  moist  soils  the  thawing 
layer  is  strong  enough  that  traction  and  motion 


Table  2.  Soil  wetness  descriptions. 


Soil  wetness 

Water  content 
(%) 

Saturation 

(%) 

wet* 

23-27 

83-100 

moist 

17-21 

77-82 

dry+ 

6-13 

23-59 

’  Nearly  saturated  and  above  the  liquid  limit. 
^  Below  the  specific  retention  of  the  soil. 


resistance  are  nearly  constant  for  all  thaw  depths 
S.  For  "dry"  soils  the  thawing  has  a  small  effect 
on  the  overall  soil  strength,  and  traction  is  a  func¬ 
tion  of  both  thaw  depth  and  applied  stress.  Since 
this  function  has  not  yet  been  defined,  we  cur¬ 
rently  use  an  interim  value  of  1.  For  "dry"  soils 
the  effect  of  thawing  on  motion  resistance  is  neg¬ 
ligible. 

The  following  expressions  define  the  multipli¬ 
ers  /and  g  in  eq  16  and  17: 

Traction: 

Wet  soils 

/=1.0  S<2.5cm  (18) 

/=  2.379  (1/S2)  +  0.619  2.5  <S<  15(19) 

/=0.63  15  <S  (20) 
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Resistance  Coefficient  Traction  Coefficient 


Thaw  Depth  (cm)  b.  Resistance  coefficient  vs.  thaw  depth. 

Figure  15.  Effect  of  soil  thaiu  depth  on  net  traction  and  motion  resistance  coefficients  as  a  function  of  soil 
moisture  content. 


Moist  soils 

/=  1.0  (21) 

Dry  soils 

/=  1.0.  (22) 

Motion  resistance: 

Wet  soils 

g  =  1.0  S<2.5cm  (23) 

y  =  2.883  S  -  6.056  2.5  <S  <4  (24) 

g  =  -0.22  S2  +  3.54  S  -  5.24  4  <  S  <  8  (25) 

g  =  0.225  S  +  7.2167  8  <  S  <  12  (26) 

g  =  10.00  12  <  S  (27) 

Moist  or  dry  soils 

g  =  1.0.  (28) 


These  equations  are  based  on  data  obtained  us¬ 
ing  an  instrumented  Jeep  Cherokee  (gross  vehicle 
weight  of  5600  lb.)  equipped  with  light-truck,  all- 
season,  radial  tires  at  a  variety  of  tire  inflation  pres¬ 
sures  (Shoop  et  al.  1993).  For  very  heavy  vehicles 
and  for  tracked  vehicles,  the  critical  thaw  depth 
will  change  because  of  the  larger  and  deeper  vol¬ 
ume  affected  by  the  applied  stress.  Because  of  the 
lack  of  effect  of  tire  inflation  pressure  (contact  pres¬ 
sure)  for  the  most  critical  conditions  (wet),  the  pro¬ 
posed  relationships  can  be  applied  to  a  wide  vari¬ 
ety  of  vehicles  with  the  understanding  that  there 
are  no  experimental  data  to  validate  the  results. 
This  aspect  will  be  addressed  in  future  studies. 

The  traction  relationships  are  based  on  the  op¬ 
timum  available  traction,  obtained  by  taking  an 
average  of  the  top  20%  of  the  traction  curve.  The 
relative  slip  values  and  traction  curve  shapes  were 
not  evaluated. 
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Because  these  relationships  were  derived  from 
tests  on  a  non-plastic  fine  sand  (unified  soil  clas¬ 
sification  SM)  (Fig.  16),  and  limited  tests  on  other 
sandy  soils  indicate  the  same  general  trends,  we 
feel  confident  applying  them  to  the  following  soil 
types: 

SW:  Well-graded  sands,  gravely  sands,  little 
or  no  fines 

SP:  Poorly  graded  sands,  gravely  sands,  little 
or  no  fines 

SM:  Silty  sands,  sand  and  silt  mixtures. 

It  is  also  reasonable  to  extrapolate  beyond  the  range 
of  test  soils  to  other  soils  of  similar  mechanical 
behavior,  such  as  other  sands  (SC),  gravels  (GW, 
GP,  GM,  GC)  and  silts  (ML).  We  don't  advise  that 
they  be  applied  to  clayey  soils  or  highly  plastic 
soils  since  the  behavior  of  the  soil,  and  there¬ 
fore  its  influence  on  traction  and  motion  resis¬ 
tance,  is  likely  to  be  quite  different.  When  thaw¬ 


ing  conditions  are  encountered  in  clays  and  me¬ 
dium-  or  high-plasticity  soils,  we  recommend  that 
algorithms  designed  to  treat  surface  slipperiness 
(Willoughby  et  al.  1991)  be  implemented  in  fu¬ 
ture  versions  of  NRMM/CAMMS.  Continuing 
research  is  aimed  at  completing  the  definition  of 
/and  g  for  more  soil  types  and  wetness  conditions. 

SPEED  MADE  GOOD 

The  time  to  get  from  one  point  to  another,  or 
speed  made  good,  is  a  function  of  many  vehicle  / 
operator  and  terrain/environment  variables.  Some 
reduction  in  speed  occurs  as  a  wheel  or  track  slips. 
Making  the  assumptions  that 

•  For  a  vehicle  moving  at  a  constant  velocity 
on  level  terrain  the  traction  generated  must 
equal  the  total  motion  resistance. 


U.S.  Standard  Sieve  Numbers  Hydrometer 


6"  4"  2"  1"  3/8-  10  40  200 


Cobbles 


Gravel 

Sand 

1  Coarse 

Fine 

C'rse 

Medium 

Fine 

Silt  or  Clay 


a.  Grain  size  distribution. 


b.  Soil  compaction  curve. 


Figure  16.  Characteristics  of 
Lebanon  sand. 
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Slip 


Figure  17.  Idealized  traction-slip  curve  in  snow. 


Figure  18.  Generalized  traction-slip  curves  for 
different  surface  materials. 


•  The  maximum  traction  in  snow  occurs  at 
about  20%  slip,  and 

•  The  traction-slip  curve  is  nearly  linear  be¬ 
tween  0  and  20%  slip  (Fig.  17), 

we  claim  that  the  slip  when  operating  in  snow 
can  be  calculated  from  the  following  equation: 


ni:..  0-2 (f^terrain  ^internal)  /•on\ 

SHPsnow  = - T - 


‘gross 


foi'  Tgross  —  ^terrain  ^internal’ 


Once  the  slip  has  been  calculated,  the  estimated 
vehicle  velocity  can  be  calculated  using  the  defini¬ 
tion  for  slip: 


Vs=Ws 


Slips 


+  1 


(30) 


where  and  Wg  are  the  vehicle  and  wheel  speeds. 


respectively.  We  recommend  that  this  velocity  cor¬ 
rection  be  applied  before  other  effects  are  taken  into 
account.  Generalized  traction-slip  curves  for  differ¬ 
ent  materials  are  shown  in  Figure  18.  These  varia¬ 
tions  win  need  to  be  considered  when  or  if  the  fuU 
traction  curves  are  incorporated  into  future  mobil¬ 
ity  models. 


INTERNAL  MOTION 
RESISTANCE 

Internal  motion  resistance,  as  mentioned  ear¬ 
lier,  is  the  force  caused  by  friction  within  the  ve¬ 
hicle  (tire  deformation,  friction  in  driveline  com¬ 
ponents,  etc.).  We  use  the  values  for  ^internal  of 
wheeled  and  tracked  vehicles  as  tabulated  by 
Ahlvin  and  Haley  (1992).  These  values  (Table  3)  are 
in  coefficient  form  (the  ratio  of  resistance  force  to 
wheel  or  track  weight). 
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Table  3.  Internal  motion  resistance  coefficients. 


^inlcrnnl 


Wheeled 

Tracked 

Surface 

0.0150 

0.0375 

Superhighways  and  primary  roads  and  ice 

0.0250 

0.0450 

Secondary  roads 

0.0175 

0.0525 

Trails  and  cross  country  (Pfg'*  >  4.0) 

0.0150 

0.0525 

Trails  and  cross  country  (Pfg  <  4.0) 

*  Pfg  =  w/ (nbr),  in  psi,  where,  w  is  the  weight  on  an  axle,  n  is  the 
number  of  wheels  on  an  axle,  b  is  the  wheel  width  and  r  is  the 
wheel  radius. 


SLOPES 

The  main  effect  of  slopes  is  the  redistribution 
of  the  force  of  gravity  on  the  vehicle  body  forces 
(Yong  et  al.  1984).  This  redistribution  is  seen  as  an 
additional  resistance  term  due  to  gravity  and  a  re¬ 
duction  of  the  normal  load  on  each  wheel  or  track 
assembly.  Richmond  (1993)  examined  this  ap¬ 
proach  using  the  traction  and  resistance  equations 
for  snow  presented  earlier  and  obtained  good 
results  when  compared  with  go /no-go  tests  on 
slopes.  The  effect  of  slopes  on  traction  is  imple¬ 
mented  by  multiplying  the  normal  load  on  a  wheel 
or  track  by  the  cosine  of  the  slope  angle.  The  effect 
on  resistance  is  implemented  by  adding  a  resistance 
term  due  to  gravity,  which  is  equal  to  the  gross  ve¬ 
hicle  weight  times  the  sine  of  the  slope  angle,  to 
^terrain' 

SUMMARY  AND 
RECOMMENDATIONS 

The  algorithms,  rationale  and  data  sources  for 
what  we  call  the  Cold  Regions  Mobility  Model 
have  been  presented.  These  algorithms  are  pri¬ 
marily  empirical  in  nature  and  rely  on  relatively 
simple  terrain  measurements.  Appendix  A  con¬ 
tains  a  stand-alone  FORTRAN  code  of  the  Cold 
Regions  Mobility  Model  (CRM-1. F);  sample  data 
files  and  output  are  presented. 

We  did  not  address  the  accuracy  of  the  algo¬ 
rithms  in  this  report;  the  reports  in  which  these 
algorithms  were  initially  developed  discuss  this 
issue  somewhat.  Fature  work  will  examine  im¬ 
provements  in  the  algorithms  and  the  overall  ac¬ 
curacy  with  regard  to  speed  made  good  and  go/ 
no-go  predictions. 


Because  this  is  an  empirically  based  model,  and 
because  mobility  is  a  continuing  area  of  interest 
to  the  U.S.  Army,  improvements  and  expansion 
of  the  model  are  to  be  expected.  In  the  process  of 
preparing  this  report  a  number  of  areas  needing 
improvement  were  identified.  For  mobility /traf- 
ficability  analysis  over  snow-covered  areas,  the  fol¬ 
lowing  need  further  research:  motion  resistance 
in  shallow  snow,  deep  snow  trafficability,  and  snow 
over  a  soft  soil  with  large  sinkage  values.  For  op¬ 
erations  on  ice  the  effects  of  temperature  and  trac¬ 
tion  aids  needs  further  study.  In  freezing/thawing 
soils  further  tests  with  heavy  vehicles  are  required, 
the  types  of  soils  investigated  need  to  be  expanded, 
and  tracked  vehicles  need  to  be  evaluated. 
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APPENDIX  A:  COLD  REGIONS  MOBILITY  MODEL  CRM-l.F 


GENERAL 


The  program  CRM-l.F  is  a  FORTRAN  code  that  calculates  mobility  parameters  and  factors  that  can  be 
used  to  determine  go/no-go  conditions,  drawbar  pull,  drawbar  coefficient,  gross  traction,  total  resistance 
and  slip  due  to  traction  on  snow.  For  calculations  on  snow  the  required  data  are  snow  depth  and  density, 
slope,  RCI  if  there  is  a  soft  soil  underlying  the  snow,  and  vehicle  data  (described  below).  If  the  analysis 
involves  only  thawing  or  freezing  soils,  the  model  requires  soil  type,  soil  wetness  (Table  2)  and  the  frost  or 
thaw  depth.  To  calculate  the  frost  depth  needed  to  support  a  vehicle,  the  vehicle  weight  or  class  is  needed. 

The  code  consists  of  a  main  program  and  the  following  subroutines: 


For  snow  calculations: 

vehdata  -  reads  the  vehicle  data  files  for  snow  modeling 

snow  -  reads  the  snow  and  terrain  data  (interactively) 

pressure  -  calculates  vehicle  contact  pressures 

density  -  calculates  pf  and  the  vehicle  sinkage 

rterrain  -  calculates  the  resistance  due  to  the  terrain 

hard  -  determines  the  internal  resistance  coefficient 

traction  -  calculates  the  gross  traction 

totalsnow  -  calculates  the  ratio  of  vehicle  to  wheel  speed,  net  traction  and  drawbar  coefficients. 


For  thawing  and  freezing  soil: 

soil  -  reads  USCS  soil  type  and  soil  wetness  (interactively) 

trmult  -  calculates  traction  multipliers  for  thawing  soils 

mrmult  -  calculates  resistance  multipliers  for  thawing  soils 

thaw  -  writes  traction  and  resistance  multipliers  to  the  screen 

freeze  -  models  frozen  soil  bearing  capacities  for  wet  or  dry  soils. 


VEHICLE  DATA  FILES 

CRM-l.F  reads  vehicle  data  files  developed  specifically  for  the  Cold  Regions  Mobility  Model  to  deter¬ 
mine  the  mobility  parameters  of  a  specified  vehicle  operating  in  snow.  No  vehicle  data  file  is  required  for 
the  freezing-thawing  soils  subroutines.  The  program  reads  two  types  of  data  files:  wheeled  and  tracked 
vehicle  files.  The  data  below  are  for  the  HMMWV,  a  wheeled  vehicle. 

The  second  line  represents  the  type  of  vehicle.  Wheeled  vehicles  are  labeled  as  type  0,  and  tracked 
vehicles  are  type  1.  Line  3  lists  the  number  of  axles.  The  number  of  wheels  per  axle  is  in  line  5;  as  indi¬ 
cated,  the  HMMWV  has  two  wheels  for  each  axle.  The  axles  are  powered  if  the  value  for  the  axle  power 
is  1;  if  the  value  is  0,  the  axle  is  not  powered.  Line  7  is  the  radius  of  each  wheel.  Line  8  is  the  maximum 


Line 

Value 

Data 

FORTRAN  format  code 

1 

Vehicle  name 

HMMWV 

2a4 

2 

Vehicle  type 

0 

10i2 

3 

No.  of  axles 

2 

i2 

4 

Ground  clearance,  inches 

15.75 

f8.2 

5 

No.  wheels  per  axle 

2 

2 

10i2 

6 

Axle  power 

1 

1 

10i2 

7 

Radius  of  each  wheel  (in.) 

16.90 

16.90 

10f7.2 

8 

Width  of  each  wheel  (in.) 

12.60 

13.00 

10f7.2 

9 

Load  on  each  axle  divided 
by  the  number  of  wheels  (lb) 

1625.00 

2135.00 

10f7.2 

10 

Tire  contact  area  (in.^) 

115.00 

115.00 

10f7.2 

11 

Unloaded  section  height 
and  loaded  tire  deflection  (in.) 

9.0 

1.2 

* 
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deflected  tire  width.  The  contact  area,  line  10,  is  the  hard  surface  contact  area  for  a  single  tire  on  each 
axle.  The  unloaded  section  height  is  determined  by  the  height  of  the  tire  without  any  load  placed  on  it  to 
cause  a  deflection.  The  loaded  tire  deflection  is  determined  by  calculating  the  deflection  of  the  tire  when 
a  load  is  placed  on  it.  Note  that  lines  8, 10  and  11  are  fimctions  of  the  tire  inflation  pressure. 

For  the  tracked  vehicles  the  data  files  are  set  up  a  little  differently.  The  example  below  is  for  the  M113. 


Line 

Value 

Data 

FORTRAN  format  code 

1 

Vehicle  name 

M113 

2a4 

2 

Vehicle  type 

1 

10i2 

3 

No.  of  track  sets 

1 

i2 

4 

Width  of  the  track  (in.) 

15.00 

10f7.2 

5 

Length  of  the  track  (in.) 

105.00 

10f7.2 

6 

Gross  vehicle  weight  (lb) 

23401.00 

f8.2 

7 

Ground  clearance  (in.) 

16.93 

f8.2 

8 

Entry  angle  (degrees) 

21.00 

10f7.2 

The  third  line  gives  the  number  of  track  sets  for  the  vehicle  (the  SUSV,  a  two-unit  articulated  vehicle  has 
two  track  sets).  The  number  of  tracks  per  set  is  set  to  2  in  the  program.  The  entry  angle,  although  re¬ 
quired,  is  not  used;  an  average  value  of  26°  for  all  tracked  vehicles  is  specified  in  CRM-1. 


TERRAIN  INPUT 

Terrain  data  are  input  interactively  through  a  series  of  menus.  For  a  soils  analysis  the  codes  in  Table 
A1  and  Table  2  (wetness  definitions)  are  used.  Table  A2  contains  the  input  data  and  associated  output  for 
a  freezing-thawing  soils  analysis.  Sample  input  data  and  the  resulting  output  are  shown  in  Table  A3  and 
A4  for  thawing  and  freezing  soils,  respectively.  For  snow  analysis  the  input  data  and  corresponding 
output  are  in  Table  A5.  Table  A6  is  a  set  of  results  for  an  analysis  of  oversnow  mobility  of  the  M113  and 
FIMMWV  (using  the  vehicle  parameters  stated  in  describing  the  vehicle  data  files),  comparing  each  of 
the  underlying  surfaces  and  snow  types. 


Table  A2.  Input  and  output  of  CRM-1  (Soils). 

Input  Output 


Thawing  Soils 


Soil  type  (Table  Al) 

Wetness  (dry,  moist,  or  wet) 
Wheeled  or  tracked  vehicle 
Thaw  depth  (cm)  _ 


yields 


Traction  multiplier 
Resistance  multiplier 


Freezing  Soils 


Table  Al.  Codes  for  different  soil  and 
terrain  surface  types  used  in  CRM-1. 
These  are  the  same  surface  codes 
used  in  CAMMS  and  NRMM. 


1 

SW 

11 

OL 

2 

SP 

12 

OH 

3 

SM 

13 

WATER 

4 

SC 

14 

PAVEMENT 

5 

SMSC 

15 

ROCK 

6 

CL 

16 

GW 

7 

ML 

17 

GP 

8 

CLML 

18 

GM 

9 

CH 

19 

GC 

10 

MH 

20 

PT  (peat) 

Wet  or  dry  ground,  plus 

Frost  depth  (cm) 

yields 

Vehicle  class 

yields 

Vehicle  weight  (lb) 

yields 

Frost  depth  and  vehicle 

yields 

Max.  vehicle  weight  and  class  supported 
Min.  frost  depth  required 
Min.  frost  depth  required 
Go/no-go 
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Table  A3.  Traction  and  motion  resistance  multipli¬ 
ers  for  wheeled  vehicles  on  thawing  sandy  soils.*  Table  A4.  Vehicle  weights  supported  by  various  frost 

depths. 


Soil  Type:  SW 

Thaw  depth  (cm) 

Wetness:  Wet 

Traction  multiplier 

Vehicle:  Wheeled 

Resistance  multiplier 

Frost  depth  (cm) 

Weight  on 
(MN) 

wet  ground 
(lb) 

Weight  on 
(MN) 

dry  ground 
'  “  (lb) 

1 

1.0 

1.0 

2 

0.0003 

77 

0.0001 

31 

2 

1.0 

1.0 

5 

0.0022 

483 

0.0008 

197 

3 

0.88 

2.6 

10 

0.0086 

1933 

0.0035 

787 

4 

0.77 

5.4 

15 

0.0194 

4350 

0.0079 

1770 

5 

0.73 

6.96 

20 

0.0344 

7734 

0.0140 

3147 

6 

0.68 

8.1 

30 

0.0774 

17401 

0.0315 

7082 

8 

0.66 

9.0 

40 

0.1376 

30935 

0.0560 

12590 

10 

0.64 

9.5 

50 

0.2150 

48336 

0.0875 

19672 

12 

0.636 

9.9 

15 

0.63 

10.0 

20 

0.63 

10.0 

*The  same  multipliers  are  obtained  (with  warnings)  for 
clayey  soils  and  tracked  vehicles.  A  multiplier  of  1.0  is 
specified  for  rock,  water  and  pavement  and  for  all  dry 
(with  warning)  or  moist  soils. 


Table  A5.  Sample  input  and  output  of  CRM-1  (snow). 


Example 

1 

2 

3 

4 

5 

Input 

Modeling  option* 

1 

1 

1 

3 

4 

'Vehicle'*' 

M113 

M113 

M113 

M113 

M113 

Snow  Base  Code** 

1 

1 

1 

3 

4 

Snow  depth  (cm) 

15 

15 

15 

15 

12 

Snow  type'*"'' 

7 

7 

6 

6 

7 

RCI 

*** 

— 

— 

— 

— 

Snow  density  (g/cm*^) 

0.2 

0.2 

— 

— 

0.2 

Slope  (degrees) 

0.0 

30.0 

0.0 

0.0 

0.0 

Output 

Vehicle  speed /wheel  speed 

0.975 

0 

0.974 

0.964 

0.94 

Rut  depth''"''''  (cm) 

9.0 

9.0 

0.0 

0.0 

7.2 

Gross  traction  (N) 

44134.3 

39207.1 

29692.1 

29692.1 

16741.3 

Total  resistance  (N) 

5585.4 

57631.6 

3903.5 

5464.8 

5275.1 

Drawbar  (N) 

38548.9 

-18424.5 

25788.6 

24227.2 

11466.2 

Drawbar  coefficient 

0.370 

-0.177 

0.248 

0.233 

0.11 

Notes 

^Modeling  option  1  is  snow  analysis,  2  is  freezing /thawing  soils. 

^Requires  a  vehicle  data  file  name. 

‘■^Snow  base  codes  (the  material  underlying  the  snow): 

1 —  superhighways /primary  roads. 

2 —  secondary  roads. 

3 —  off-road  terrain/trails/strong  frozen  ground. 

4 —  ice  (if  snow  depth  is  <1.0  cm,  then  ice  only). 

5 —  unfrozen  soil  (requires  RCI). 

■•^'''Snow  type: 

6 —  disturbed  (more  than  four  vehicle  passes  or  a  prepared  snow  road). 

7 —  undisturbed  snow. 

***Input  not  required  with  the  previous  specified  input. 

'’"’"''Total  vehicle  sinkage. 
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Table  A6.  Terrain/vehicle  data  sets  from  CRM-1. 


Vehicle:  M113 


IBASE 

Vs/Ws 

Rut  depth  (cm) 

Gross 

traction  (N) 

Total 

resistance  (N) 

Drazubar  (N) 

Drawbar 

coefficient 

Snow  Type: 
Snow  Depth: 

Undisturbed 
15  cm 

Density: 

0.2  Mglm^ 

Slope: 

0,0° 

1 

0.975 

9.0 

44134.3 

5585.5 

38548.9 

0.370 

2 

0,972 

9.0 

44134.3 

6366.1 

37768.2 

0.362 

3 

0.968 

9.0 

44134.3 

7146,8 

36987.5 

0.355 

4 

0.937 

9.0 

16741.3 

5585,5 

11155.9 

0.107 

5  (RCI=30) 

0.965 

14.1 

44134.3 

7996.6 

36137.7 

0.347 

Snow  Depth: 

20  cm 

Density: 

0.15  Mglm^ 

Slope: 

10.0° 

1 

0.901 

14.0 

43581.8 

23915.4 

19666.4 

0.188 

Snow  Depth: 

0.0 

Density: 

— 

Slope: 

0.0° 

4 

0.930 

0.0 

10409.2 

3903.5 

6505.8 

0.063 

Snow  Type: 
Snow  Depth: 

Disturbed 

15  cm 

Density: 

— 

Slope: 

0.0° 

1 

0.974 

0.0 

29692.1 

3903.5 

25788.6 

0.248. 

2 

0.969 

0.0 

29692.1 

4684.2 

25007.9 

0.240 

3 

0.964 

0.0 

29692.1 

5464.8 

24227.2 

0.234 

4 

0.974 

0.0 

29692.1 

3903,5 

25788.6 

0.248 

5  (RCI=20) 

0.935 

10.5 

29692.1 

10359.0 

19333.0 

0.186 

Snow  Depth: 

15  cm 

Density: 

— 

Slope: 

15.0° 

1 

0.0 

0.0 

28710.2 

30844.5 

-2134.4 

-0.021 

Snow  type: 
Snow  depth: 

Undisturbed 
15  cm 

Density: 

0.2  Mg/m^ 

Slope: 

0.0° 

1 

0.934 

9.0 

12318.2 

4348,9 

7969.3 

0.238 

2 

0.929 

9.0 

12318.2 

4683.4 

7634.8 

0.228 

3 

0.933 

9.0 

12318.2 

4432.5 

7885.7 

0.236 

4 

0.866 

9.0 

5643.9 

4348.9 

1295.0 

0.039 

5  (RCI=30) 

0.924 

12.6 

12318.2 

5096.3 

7221.9 

0.216 

Snow  Depth: 

20  cm 

Density: 

0.15  Mglm^ 

Slope: 

10.0° 

1 

0.859 

14.0 

12163.9 

9965.6 

2198.4 

0.068 

Snow  Depth: 

0.0 

Density: 

— 

Slope: 

0,0° 

4 

0.971 

0.0 

3345.0 

501,8 

2843.3 

0.085 

Snow  Type: 
Snow  Depth: 

Disturbed 

15  cm 

Density: 

— 

Slope: 

0.0° 

1 

0.989 

0.0 

9316.0 

501.8 

8814.3 

0.264. 

2 

0.982 

0,0 

9316.0 

836.3 

8479.8 

0.254 

3 

0.986 

0.0 

9316.0 

858.4 

8730.7 

0.261 

4 

0.989 

0.0 

9316.0 

501.8 

8814.3 

0.264 

5  (RCI=20) 

0.835 

7.0 

9316.0 

9193.4 

122.6 

0.004 

Snow  Depth: 

15  cm 

Density: 

— 

Slope: 

15.0° 

1 

0.0 

0.0 

9008.0 

9159.4 

-151.4 

-0.005 
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C  CRM-l.f 
C  7  MARCH  1994 


C  THIS  PROGRAM  CALCULATES  THE  TRACTION  AND  MOTION  RESISTANCE 

C  OF  WHEELED  AND  TRACKED  VEHICLES  OPERATING  ON  ICE,  OR  IN  DEEP  OR 

C  SHALLOW  SNOW  OVER  VARIOUS  BASE  MATERIALS  AND  THAWING  SOILS 


C  PROGRAM  BY: 

C  HEATHER  BANKER 

C  GEORGE  BLAISDELL 

C  PAUL  RICHMOND 

C  SALLY  SHOOP 

C  KAREN  FARAN 

common/ml/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 

common/m5/iveh,nunits,pmax,gclear,gvw 

common/m6/ def,sec 

common/m7/ibase,rci,idist,slope 

1000  write(*,*)'choose  an  option:' 

write(*,*)'l  =  snow  modeling  (snow  or  ice  covered)' 
write(*,*)'2  =  soils  modeling  (frozen  or  thawing  soils)' 
read(*,*)input 

C  INITIALIZE  VARIABLES 
sumz=0. 
zmax=0. 
gvw=0. 
rsnow=0.0 
trsnow=0. 
totalr=0. 
totaldb=0. 

if(input.eq.l)then 
call  vehdata 

call  snow(hsnow,sigmaO) 
call  pressure 

call  density(hsnow,sigmaO,sigmaf,zmax) 

call  rterrain(hsnow,sigmaO,sigmaf,zmax,sumz,rsnow) 

call  hard(reshard) 

call  traction(hsnow,trsnow) 

call  totalsnow(rsnow,trsnow,reshard,hsnow,sigmaO,sumz) 
else 

write(*,*)'what  type  of  soil  problem  is  this' 
write(*,*)'l  =  thawing  soil' 
write(*,*)'2  =  frozen  soil' 
read(*,*)iput 

if(iput.eq.l)then 
call  soil(isoil,iwet) 

call  trmult(isoil,iwet,trmlt,zthaw,ivehs) 

call  mrtmult(zthaw,svalid,isoil,iwet,rmult,ivehs) 
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call  thaw(trmlt,rmult,zsnow) 
else 

call  freeze(isoil) 
end  if 
end  if 

write  (’^/)'enter  y  to  run  another  condition' 
read(%300)ans 
300  format  (al) 

if  ((ans.eq."Y")  .or.  (ans.eq."y"))  then 
goto  1000 
end  if 

c  call  exit(O) 

end 


subroutine  vehdata 


C  THIS  SUBROUTINE  READS  DATA  FROM  THE  VEHICLE  DATA  FILE 

common/ml/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 

common  /  m5  /  iveh,nunits,pmax,gclear,gv  w 

common/  m6  /  def  ,sec 

character’^12  mfile 


write(*,*)'enter  the  vehicle  data  file  name' 
read(*,10)  mfile 

open(21,status='old',file=mfile,form='formatted') 

C  READ  IN  VEHICLE  DATA 
read(21,100)title 

read(21,110)iveh  !WHEELED=0,  TRACKED=1 


if(iveh.eq.O)  then 


C  READ  WHEEL  VEHICLE  DATA  • 
read(21,115)nunits 
read(21,220)gclear 
read  (21, 11 0)  (nwhls  (i),i= l,nunits) 
read  (21, 110)  (ip  (i),i= l,nunits) 
read(21,120)(rads(i),i=l,nunits) 
read(21,120)(width(i),i=l,nunits) 
read(21,120)(pload(i),i=l,nunits) 
read(21,120)(area(i),i=l,nunits) 
read(21,‘^)sec,def 
TIRE  DEFLECTION 


ENGLISH  UNITS 
!NO.  AXLES, 

IGROUND  CLEARANCE 
INUMBER  OF  WHEELS  PER  AXLE 
lAXLE,  0=N0  POWER  ,l=POWERED 
IRADIUS  OF  EACH  WHEEL,  IN. 

(WIDTH  OF  EACH  WHEEL,  IN. 

ILOAD  ON  EACH  WHEEL,  LBS. 

ITIRE  CONTACT  AREA,  SQ.  IN. 
[UNLOADED  SECTION  HEIGHT,  LOADED 
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C  CALCULATE  GVW  FOR  WHEELED  VEHICLES 
do  i=l,nunits 

gvw=gvw+(pload(i)*nwhls(i)) 
end  do 

write(L115)nunits 

write(*,220)gclear 

write(*,110)(nwhls(i),i=l,nunits) 

write(L110)(ip(i),i=l/nunits) 

write(L120)(rads(i),i=l,nunits) 

write(L120)(width(i),i=l,nunits) 

write(L120)(pload(i),i=l,nunits) 

write(L120)(area(i),i=l,nunits) 

write(L’^)sec,def 

write  (*/)gvw 

else 

C  READ  TRACK  VEHICLE  DATA  -  ENGLISH  UNITS 
read(21,215)nunits 

read(21,120)(width(i),i=l/nunits)  ITRACK  WIDTH 

read(21,120)(tlength(i),i=l,nunits)  ITRACK  LENGTH 

read(21,220)gvw  IGROSS  VEHICLE  WEIGHT 

read(21,220)gclear  IGROUND  CLEARANCE 

read(21,120)(aa(i),i=l,nunits)  lENTRY  ANGLE 

write(*,215)nunits 

write(*,120)(width(i),i=l,nunits) 

write(L120)(tlength(i),i=l,nunits) 

write(*,220)gvw 

write(L220)gclear 

write(L120)(aa(i),i=Lnunits) 

end  if 

10  format(al2) 

100  format(2a4) 

110  format(10i2) 

115  format(i2) 

120  format(10f7.2) 

130  format(f6.2,f6.3) 

215  format(i2) 

220  format(fl0.2) 

close(21) 

return 

end 

subroutine  snow(hsnow,sigma0) 

C  THIS  SUBROUTINE  READS  IN  SNOW  AND  TERRAIN  DATA 
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common  /  m7/ ibase,rci,idist,slope 

C  ENTER  TERRAIN  DATA 

write(^*)'enter  snow  base  code:' 
write(^*)'l  =  super  highways/primary  roads' 
write(*,*)'2  =  secondary  roads' 

write('*^/)'3  =  off-road  terrain/trails/strong  frozen  soil' 
write(*/)'4  =  ice  (if  snow  depth  is  <  1.0cm  then  ice  only)' 
write(*,*)'5  =  unfrozen  soil' 
read(*/)  ibase 

write(’'^/)'input  the  snow  depth  (cm)' 
read(*,*)hsnow 

if  ((ibase.ne.4).or.(hsnow  .ge.  1.0))  then 
write(*,*)'enter  the  type  of  snow' 

write(’^/)'6  =  disturbed  snow  (more  than  4  vehicle  passes)' 

write(*,*)'7  =  undisturbed  snow' 

read(*/)idist 

end  if 

if  (ibase.eq.5)  then 

write(*,*)'mput  rating  cone  index  of  the  soil' 
writer, *)'20  <  RCI  <  100' 
read(*,*)  rci 
end  if 

if  (idist.eq.6)  then 
if  (sigmaO.le.0.55)  then 
sigmaO  =  0.55 
end  if 
else 

if  ((ibase. ne.  4).or.(hsnow  .ge.  1.0))  then 
write(*,*)'mput  the  snow  density  (g/cm**3)' 
read('^,*)sigma0 

end  if 
end  if 

write(^'^)'enter  the  angle  of  the  slope  (degrees)' 
read(*/)slope 

return 

end 


subroutine  pressure 

C  THIS  SUBROUTINE  CALCULATES  THE  CONTACT  PRESSURE 

common/ml/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 
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common  /  m5  /  iveh,nunits,pmax,gclear,gvw 
if  (iveh.eq.l)  then 

C  CALCULATE  VALUES  EOR  TRACKED  VEHICLES 
do  i=l,nunits 

nwhls(i)=2  !  ALL  TRACKED  VEHICLES  HAVE  2  TRACKS  PER  AXLE 
ip(i)=l  !ALL  TRACKS  ARE  DRIVEN 

area(i)=width(i)*tlength(i) 
pIoad(i)=gvw/nunits/nwhIs(i) 
presl(i)=pload(i)/ area(i) 
c  write(L*)'pressure',presl(i) 

end  do 

else 

C  CALCULATE  GROUND  CONTACT  PRESSURES  FOR  WHEELED  VEHICLES 
do  i=l,nunits 
presl(i)=pload(i)/area(i) 
c  write(L*)'pressurel',presl(i) 

end  do 
end  if 

C  CALCULATE  THE  MAXIMUM  PRESSURE 

pmax=presl(l) 
do  i=l,nunits 

if  (presl(i).gt.pmax)  pmax=presl(i) 
end  do 

return 

end 

subroutine  density(hsnow,sigmaO,sigmaf,zmax) 

C  THIS  SUBROUTINE  CALCULATES  THE  FINAL  DENSITY  AND  SINKAGE 

common  /  m5  /  iveh,nunits,pmax,gclear,g  vw 

C  COMPUTE  FINAL  DENSITY  DATA 
sigmaf=.5 

if(pmax.gt.30.5)  sigmaf=.55 
if(pmax.gt.50.8)  sigmaf=.60 
if(pmax.gt. 101.5)  sigmaf=.65 

C  COMPUTE  SINKAGE 

zmax=  hsnow  *  (1.  -  (sigmaO/sigmaf)) 
if(zmax.lt.0.)zmax=0. 

return 

end 

subroutine  rterrain(hsnow,sigmaO,sigmaf,zmax,sumz,rsnow) 
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C  THIS  SUBROUTINE  CALCULATES  TERRAIN  RESISTANCE  OF  THE  VEHCILE 


common/ml /nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 

common  /  m5  /  iveh,nunits,pmax,gclear,gvw 

common /m6/def, sec 

common/ m7/ibase,rci,idist,slope 

dimension  xl2(10),rresl(10) 

C  START  LOOP  FOR  RESISTANCE 
flag  =  1.0 
rsnow  =  0.0 
soilz  =  0.0 

doi=l,10 
xI2(i)=0.0 
rresl(i)=0.0 
end  do 

i=l  !  ALL  RESISTANCE  CALCULATIONS  ARE  BASED  ON  FIRST  ! 

UNIT  CHARACTERISTICS  AND  PMAX 

C  WHEELED  VEHICLES 
if  (iveh.eq.O)  then 

C  ACCT  FOR  WEAK  SOIL  BASE 
if(ibase.eq.5)  then 

pn  =  1.  !IS  THE  NUMBER  OF  PASSES  THE  VEHICLE  MAKES 

s  =  0.05  !SLIP  ASSUMMED  CONSTANT  FOR  THIS  CALC 

C  def  IS  THE  LOADED  TIRE  DEFLECTION 
C  sec  IS  THE  UNLOADED  TIRE  SECTION  HEIGHT 

soilz=(10.'^rads(i)*pn*'’^.5)/((rci*2.*rads(i)*width(i))/ 
&(pIoad(i)^(l.-def/sec)^’^(3./2.)*s’'*0.2))**(5./3.) 
if  (soilz. It. 0.0)  soilz=0.0 

soilz  =  soiIz*2.54 

end  if 

sumz=zmax+soilz  !sumz  IS  IN  cm 

C  RESISTANCE  PARAMETER,  WHEELED  VEHICLES  IN  SNOW 

xI2(i)=rads(i)*acos((rads(i)-(sumz/2.54))/rads(i))  IINCHES 
xI2(i)=xI2(i)=^2.54  !cm 
c  write(*,'^)'xl2',xI2,'rads',rads(i) 

C  ACCOUNT  FOR  TIRE  CURVATURE/DEEP  SNOW  FLAG 
radscm  =  rads(i)'^2.54 
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if((sigma0.1t.0.15).and.(sumz.ge.radscm))flag  =  1.5 
if((sigma0.ge.0.15).and.(sumz.ge.(radscm*2./3.)))flag  =  1.5 
if((sigma0.ge.0.15).and.(suinz.ge.(radscm)))flag  =  2.5 

C  GROUND  CLEARANCE  FLAG 

if  (sumz.ge.  (gclear*2.54))  flag  =  flag+1.5 

else 

C  TRACKED  VEHICLE  RESISTANCE 
C  ACCOUNT  FOR  WEAK  SOIL  BASE 

if(ibase.eq.5)  then 

pn  =  1.  !IS  THE  NUMBER  OF  PASSES  THE  VEHICLE  MAKES 

soilz=tlength(i)*pn**.5*0. 00443* 

&exp((5.889*pload(i))/(rci*width(i)*tlength(i))) 
soilz  =  soilz*2.54 
end  if 

sumz=zmax+soilz  !cm 

C  RESISTANCE  PARAMETER,  TRACKED  VEHICLES  IN  SNOW 
aa(i)=26.0  !AN  AVERAGE  VALUE  FOR  NRMM 
xl2(i)=(sumz/sm(aa(i)*3.14159265/180.))  !cm 

C  GROUND  CLEARANCE,  LOW  (FLAG=2.0)  OR  HIGH  (FLAG=4.0)  DENSITY  SNOW 

if((sigma0.1t.0.15).and.(sumz.ge.(gclear*2.54)))  flag  =  2.5 
if((sigma0.ge.0. 15). and. (sumz.ge. (gclear*2.54)))  flag  =  4.0 

end  if 

if  (flag.gt.1.0)  then 
C  DEEP  SNOW  WARNING 
write(*,*) 

write^,*)'A  DEEP  SNOW  FLAG  HAS  BEEN  SET,' 
write(*,*)'CALCULATIONS  ARE  NOT  VALIDATED.' 
write(*,*)'FLAG  =  ',  flag 
write(*,*) 

end  if 

c  print  *,  rads(i) 

print  *,  xl2(i),sumz 
xl=xl2(i) 

w=width(i)*2.54  (CONVERT  TO  cm 

rresl(i)=(68.1*(0.10*sigma0*xl*w)**0.914)*flag  IKG/M  IN  PARATHESISES 
c  write(*,*)'flag',flag,'rresl  ',rresl,'sigmaO',sigmaO 

c  write(*,*)(0.10*sigma0*xl*w),xl,w,sigma0 
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rresl(i)=rresl(i)/4.4482  ICONVERT  TO  lbs 


C  ADD  SLOPE  TO  RESISTANCE  AND  SUM  RESISTANCE  FOR  WHOLE  VEHICLE 
C  UP  TO  4  UNITS  ONLY  FOR  WHEELED  VEHICLE  SNOW  RESISTANCE, 

C  1  UNIT  LIMIT  FOR  TRACKS 

sloper=slope*3.14159265/180.  [SLOPE  IN  RADIANS 

if(iveh.eq.O)then 
do  j=l,nunits 
rresl(j)=rresl(l) 

if  (j.gf-4)  then 
rresl(j)=0.0 
end  if 

rsnow=  rsnow+((rresl(j)+pload(j)“'sin(sloper))*nwhls(j)) 

end  do 

else 

rsnow=rresl(l)*nwhls(l)+gvw’^sin(sloper) 
end  if 

return 

end 


subroutine  hard(reshard) 

common/ml/nwhls(10),rads(10),width(10),pload{10),area(10),aa(10) 
common/m2/ip(10),presl(10),tlength(10) 
common/m5/iveh,nunits,pmax,gclear,gvw 
common  /  m6  /  def  ,sec 
common/  m7  /  ibase,rci,idist,slope 

dimension  cpffg(lO) 


C  CPFFG  =  F-G  CONTACT  PRESSURE  FACTOR  FOR  EACH  ELEMENT 
C  W  =  WEIGHT  BENEATH  SUSPENSION  ASSEMBLY  (TRACK  PAIR  OR  AXLE)  IN  LBS 
C  RHARD  =  HARD  SURFACE  RESISTANCE  COEFFICIENT 


reshard  =  0.0 

C  WHEELED  VEHICLES 
if(iveh.eq.O)then 

do  i  =  l,nunits 

w  =  pload(i)*nwhls(i) 

cpffg(i)  =  w/(nwhls(i)’^width(i)*rads(i)) 


if((ibase.eq.3).or.(ibase.eq.5))then 

if(cpffg(i).ge.4.0)then 
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rhard  =  0.0175 
else 

rhard  =  0.015 
end  if 

else  if((ibase.eq.l).or.(ibase.eq.4))then 

rhard  =  0.015 

else 

rhard  =  0.025 
end  if 

reshard  =  reshard  +rhard*w 
end  do 

C  TRACKED  VEHICLES 
else 

if((ibase.eq.3).or.(ibase.eq.5))then 
rhard  =  0.0525 

else  if((ibase.eq.l).or.(ibase.eq.4))then 
rhard  =  0.0375 

else 

rhard  =  0.045 
end  if 

reshard  =  rhard*gvw  ITOTAL  HARD  SURFACE  RESISTANCE  FOR  TRACK 
end  if 

return 

end 


subroutine  traction(hsnowTrsnow) 

C  THIS  ROUTINE  CALCULATES  THE  TRACTION  OF  THE  VEHICLE 


common/ml/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 

common  /  m5  /  iveh,nunits,pmax,gclear,gvw 

common/ m6/def,sec 

common/m7/ibase,rcLidist,slope 

dimension  trcnl(lO) 

C  TRACTION  CALCULATION 
trsnow  =  0.0 
do  i  =  l,nunits 

pr=presl(i)’^6.8948  !  KPA 

pr=pr*cos(slope^3.14159265/180.)  !  ACCOUNT  FOR  SLOPE 

C  UNDISTURBED  SNOW 

if(((ibase.eq.l).or.(ibase.eq.2).or.(ibase.eq.3).or.(ibase.eq.5)) 
&.and.(idist.eq.7))then 

trcnl(i)=0.851*((pr)=^=^0.823) 
write(L*)'area',area(i)/load',pIoad(i) 
write(L*)'pressure',presl(i) 


c 

c 
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c 

c 


write(’^/)'pressure  in  kPa',pr 
write  (’'■/)trcnl(i) 
end  if 


C  HARD  PACKED  SNOW 
if  (idist.eq.6)  then 

if((ibase.eq.l).or.(ibase.eq.2).or.(ibase.eq.3))then 

trcnl(i)=0.32r(pr)**0.97 

end  if 

if(ibase.eq.5)then 
trcnl(i)=0.321='(pr)**0.97 
end  if 
end  if 

C  SNOW  ON  ICE  BASE 

if((ibase.eq.4).and.(idist.eq.7))then 

trcnl(i)=0.127*(pr)*n.06 

else  if((ibase.eq.4).and.(idist.eq.6))then 

trcnl(i)=0.321=^(pr)*=^0.97 

end  if 

C  ICE  SURFACE 

if  ((ibase.eq.4).and.(hsnow.le.l.0))then 

trcnl(i)=O.Hpr 

end  if 

trcnl(i)=(trcnl(i)/6.8948)»area(i)  ICONVERT  TO  lbs 

C  NON-DRIVEN  WHEEL 

if(ip(i).eq.O)  trcnl(i)=0.0 

C  SUM  TRACTION  FOR  EACH  ELEMENT 
trsnow  =  trsnow  +  trcnl(i)*nwhls(i) 

end  do 

return 

end 

subroutine  totalsnow(rsnow,trsnow,reshard,hsnow,sigmaO,sumz) 

common/ml /nwhls(10),rads(10),width(10),pload(10),area(10),aa(10) 

common/m2/ip(10),presl(10),tlength(10) 

common/m3/xl2(10),cpffg(10) 

common  /  m5  /  iveh,nunits,pmax,gclear,gvw 

common/m6/ def,sec 

common  /  m7/ibase,rci,idist,slope 

C  TOTAL  VEHICLE  TRACTION  AND  RESISTANCE 
totalr  =  reshard  +  rsnow 
totaldb  =  trsnow  -  totalr 
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C  SLIP  DUE  TO  TRACTION 

if  (totalr.Ie.trsnow)  then 
slipt  =  .2*totaIr  /  trsnow 

slipw  =  1. /(slipt  +  1.)  '.RATIO  OF  VEHICLE  SPEED  TO  WHEEL  SPEED 
else 

slipw=0.0 
end  if 

C  WRITE  TOTAL  VEHICLE  TRACTION  AND  RESISTANCE 
write(%210)hsnow,sigma0,slope 

write(’^,*)'ratio  of  vehicle  speed  to  wheel  speed:  ',slipw 
write(L*)'  units  are  lbs:' 

write(L260)sumz,trsnow,totaIr,totaldb 
write(*,*)'  units  are  N:' 

write(’^,260)sumz,trsnow*4.4482,totalr*4.4482,totaldb*4.4482 
write  (*,*)'drawbar  coefficient  =  ;(totaIdb/gvw) 

200  format(lx,al2) 

210  format(/lx,"snow  depth,  cm.=  ",f5.2,2x/'density,  g/cc=  ",f5.4, 

&  "  slopes  ",f6.2) 

260  format(lx/'rut,cm.=  ",f4.1,2x,"tr  =",f8.1,2x/'rr=  ",f8.1, 

&  2x,"db=  ",f8.1) 

return 

end 

subroutine  thaw(trmIt,rmult,zsnow) 

C  if(layer.eq.l.and.zsnow.It.l.O/ xcm2in)then 

do  i  =  1,  nunits 

c  trmltn  =  (dowpb(i)  +  rtowpb(i))*trmIt 

c  rtowpb(i)  =  rtowpb(i)*rmult 

c  rtowt(i)  =  rtowt(i)*rmuIt 

c  dowpb(i)  =  trmltn  -  rtowpb(i) 

end  do 
C  else 

C  write(’^,*)'routine  thaw  is  not  implemented  for:' 

C  if(layer.ne.l.)then 

C  write(L*)'critical  depth  >  6  in.' 

C  else 

C  write(L*)'snow  cover  >  1  cm' 

C  end  if 

C  end  if 

write(*,*)  'traction  multiplier  ',trmlt 
write(*/)'motion  resistance  multiplier ',rmult 
C 

return 

end 


subroutine  trmult(isoil,iwet,trmlt,zthaw,ivehs) 
c  this  subroutine  calculates  the  traction  multiplier 
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write(%*)'are  you  using  wheeled  (0)  or  a  tracked  vehicle  (1)' 
read(’^,*)ivehs 


write(^*)'input  the  thaw  depth  (cm)' 
read(’^/)zthaw 

trmlt=  1. 

if((iwet.eq.3).and.((isoil.lt.l3).or.(isoil.gt.l5)))then 
if(zthaw.le.2.5)then 
trmlt  =  trmlt 

else  if((zthaw.gt.2.5).and.(zthaw.lt.l5.))then 
trmlt  =  trmlh^(2.379’^(l/zthaw**2.)+0.619) 
else  if(zthaw.ge.l5.)then 
trmlt  =  0.63"^trmlt 
end  if 

else  if(iwet.eq.2)  then 
trmlt  =  1. 

else  if(iwet.eq.l)then 
write 

write(*/)'for  dry  soils,  the  effect  of  thawing  on  traction  is  a' 
writeC^,*)'function  of  applied  stress;  these  functions  have  not 
write(^*)yet  been  developed  so  a  multiplier  of  1.0  is  used' 
write(*  ’*'y^*******’*'’f”*'**^*=^***’^*’*'*’^********’‘'’*'**=<-=*'*’f-’t-=t-=<-*******=^’*-=*-+*=<-*' 
end  if 

if((isoil.ge.l3).and.(isoil.le.l5))then 

trmlt  =  1. 

write 

write(*,*)'A  multiplier  of  1.0  is  used  for  water,  rock,  or  pavement' 
write  (* 
end  if 

if(((isoil.ge.6).and.(isoil.le.l2)).or.(isoil.eq.20))then 

write(*  *y^**’*'**’^*^’^’*'’*'******=*'=f'’^’^’*'’f'*’*'*****=*'=*'**’f-*=+*=*- <■=♦■*’*■’*■*=«■********' 

write(*,’^)'calculations  are  based  on  sandy  soils' 

write (* 

end  if 

if(ivehs.eq.l)then 

write (*  *y’*'’*'^***’*'’*'’^’f'*’f^’*'=*^^’*'**=t'**’*'*’i^*<'*=*-*+*’i-=*'**=t^=h*=<-**=<'=*-' 
write(*,’^)'calculations  are  based  on  wheeled  vehicles  ' 
write 
end  if 

return 

end 


subroutine  mrtmult(zthaw,svalid,isoil,iwet,rmult,ivehs) 
c  this  subroutine  calculates  the  motion  resistance  multipler 
rmult  =  1. 

c  if(iwet.eq.3)then 

if((iwet.eq.3).and.((isoil.lt.l3).or.(isoil.gt.l5)))then 
if(zthaw.le.2.5)then 
rmult  =  1. 

else  if((zthaw.gt.2.5).and.(zthaw.lt.4.))then 
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rmult  =  rmult’*'(2.883*zthaw  -  6.056) 
else  if((zthaw.ge.4.).and.(zthaw.le.8.))then 

rmult  =  rmult*(-0.22*zthaw**2.  +  3.54*zthaw  -  5.24) 
else  if((zthaw.gt.8.).and.(zthaw.lt.l2.))then 
rmult  =  rmult'^(0.225*zthaw  +  7.2167) 
else  if(zthaw.ge.l2)then 
rmult  =  rmulO'lO.O 
end  if 
end  if 

return 

end 


subroutine  soil(isoil,iwet) 

c  this  subroutine  asks  the  USCS  soil  type  and  wetness 
write(*/)'choose  an  USCS  soil  type  code:' 
write(*/)'  1  =  SW' 
write(^*)'  2  =  SP' 
write(*/)'  3  =  SM' 
write(*/)'  4  =  SC' 
write(*/)'  5  =  SMSC' 
write(^*)'  6  =  CL' 
write(*,»)'  7  =  ML' 
write(*,*)'  8  =  CLML' 
write(*,*)'  9  =  CH' 
write(*/)'10  =  MH' 
write(*/)'ll  =  OL' 
write(*/)'12  =  OH' 
write(*/)'13  =  WATER' 
write(*/)'14  =  PAVEMENT' 
write(*/)'15  =  ROCK' 
write(^*)'16  =  GW' 
write(^*)'17  =  GP' 
write(*A)'18  =  GM' 
writeC/T19  =  GC' 
write(L*)'20  =  PT' 
read(*A)isoil 

write(’^A)'input  the  soil  wetness' 
write^A)'  1  =  dry' 
write(%*)'  2  =  moist' 

write(*A)'  3  =  wet  (water  content  is  above  liquid  limit)' 
read(*,*)iwet 

return 

end 

subroutine  freeze 

writeC<'A)'FROZEN  GROUND  MENU' 
write(%*)'l  =  input  the  frost  depth' 
write(*A)'2  =  input  the  vehicle  class' 
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write(*,*)'3  =  input  the  vehicle  weight  (lb.)' 

write(*/)'4  =  input  frost  depth  and  vehicle  class  or  weight' 

read(*/)nput 

write(%*)'choose  a  soil  condition' 
write(*,*)'l  =  wet' 
write(*,*)'2  =  dry' 
read(*/)cond 

if(nput.eq.l)then 

write(^*)'input  the  frost  depth  (cm)' 
read(*/)fdepth 

fdepth  -  fdepth/ 100.  ICONVERT  TO  METERS 

C  CALCULATE  THE  MAXIMUM  LOAD  THAT  A  SPECIFIED  FROST  DEPTH  WILL  HOLD 
if(cond.eq.l)then 
p  =  0.86*fdepth**2. 
else  if(cond.eq.2)then 
p  =  0.35*fdepth**2. 
end  if 

write(*,*)'maximum  load  that  can  travel  on'ddepth/m  =',p/MN' 
c  convert  dimensions  of  p  from  MN  to  pounds 
p=p/(4.448222e-6) 

write(*,*)'maximum  load  that  can  travel  on',fdepth/m  =',p/pounds' 
else  if(nput.eq.2)then 
write(*,*)'input  the  vehicle  class' 
read(*,*)vclass 

C  CALCULATE  THE  MINIMUM  FROST  DEPTH  THAT  WILL  HOLD  A  SPECIFIED  VEHICLE 
if(cond.eq.l)then 
fdepth  =  10*(vclass)**.5 
else  if(cond.eq.2)then 
fdepth  =  16*(vclass)**.5 
end  if 

fdepth  =  fdepth/ 100  (convert  from  cm  to  m 

write(*,*)'mrnimum  frost  depth  that  can  hold  class',vclass/=', 

&  fdepth,'m' 

else  if(nput.eq.3)then 

c  write(*,*)'input  the  vehicle  weight  (MN)' 

write(*,*)'input  the  vehicle  weight  (pounds)' 
read(*,*)vwght 

vwght=vwght*4.448222e-6  (convert  to  MN 

C  CALCULATE  THE  MINIMUM  FROST  DEPTH  THAT  WILL  HOLD  A  SPECIFIED  LOAD 

if  (cond  .eq.  1  )then 
c  fdepth  =  (0.86/p)=^L5 

fdepth  =  (vwght/0.86)**.5 
else  if(cond.eq.2)then 
c  fdepth  =  (0.35/p)=^L5 

fdepth  =  (vwght/0.35)*L5 
end  if 

write(*/)'mmimum  frost  depth  that  can  hoId',vwght/MN  =',fdepth/m' 
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vwght=vwght/4.448222e-6  Iconvert  to  pounds 

write(*,*)'minimum  frost  depth  that  can  hold',vwght/pounds  =', 

&  fdepth/m' 

C  CALCULATE  THE  GO/NO  GO  SITUATION  OF  BOTH  VALUES  INPUTTED 
C  EITHER  THE  WEIGHT  OR  THE  CLASS  WILL  BE  READ 


else  if(nput.eq.4)then 

write(*,*)'input  the  frost  depth  (cm)' 
read(*,*)fdepth 

fdepth  =  fdepth/100.  ICONVERT  TO  METERS 
write(*,*)'are  you  using  vehicle  weight  (1)  or  class  (2)' 
read(*,*)ans 
if(ans.eq.l)then 

write(*/)'input  the  vehicle  weight  (MN)' 
write(*,*)'input  the  vehicle  weight  in  pounds' 
read(*,*)vwght 
vwght=vwght‘^4.448222e-6 
if(cond.eq.l)then 
if(vwght.le.0.86*fdepth**2.)then 

write(*/)'the  ground  can  support  this  vehicle!' 
writ6(*  ’f'y****’*'**’^*****’*’*’*'**’*'’^*******’*'^********’*'’^’*'' 

else 

write(*,*)'the  ground  can  not  support  this  vehicle!' 
writ6(*  ’‘■y’^**’*'****’^'*’*'*****’*'’*'*********’*'’^’*”*'****’*’*’*'’*”*'' 

end  if 

else  if(cond.eq.2)then 
if(vwght.le.0.35*fdepth**2.)then 

write(’*'  ’‘■y*************’*”*'****’*'**’*'’*”*”*'*******’*'*****’*’' 
write(*/)'the  ground  can  support  this  vehicle!' 
write(*  ’(■y 

else 

write(*  ’f'y’^*^^**’*'***’''*’*'’''’*'*****’*"*'*’*'’*'*******’*'***'*'***’*'' 
write(*/)'the  ground  can  not  support  this  vehicle!' 
write(*  >f-)'***)f-*>t-sf-****=('**>ntH-*=f*****<'>t-’f'*+***=<-<-*+’t'*’t*' 

end  if 
end  if 

else  if(ans.eq.2)then 
write(*,*)'input  the  vehicle  class  (tons)' 
read(*,*)vclass 
if(cond.eq.l)then 
if  (fdepth.  ge.  10.  *(vclass)*’^.5)then 

write(’*'  ’‘■y’f'*’*'’*'’*'’^'*****’*'*’*'***’*'***'^’*'*’^**’*'’*'*’*"'”^’^*****’*'' 
write(*/)'the  ground  can  support  this  vehicle!' 
write(’*'  ^y^’f^’f"’*"’*'*’*'’^**’*'’*'’*'*’*'^*’*^*’*'**’^’*'****’*'****’*'’*'^****' 

else 

write(’*'  ^y’*”*'’^’*'**’^'*’^'’^’*'****’*'*’*'’*”^*******’*'’^**’^’^****’*'’*'*' 
write(*/)'the  ground  can  not  support  this  vehicle!' 
write(*  ’♦■y**’^'*’^'’*'’*'****’^'’*'’*'*’*'^****’*'’*'’*'*’*'’*'**’*'******’*'*’*'*' 

end  if 

else  if(cond.eq.2)then 
if(fdepth.ge.l6.*(vclass)**.5)then 
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else 


writ6(^ 

write(*/)'the  ground  can  support  this  vehicle!' 
write (* 

write (* 

write(’^/)'the  ground  can  not  support  this  vehicle!' 


end  if 
end  if 
end  if 
end  if 

if(fdepth.ge.0.5)then 
write^/)' 


5f•5f5f5f■:fsf•:f5f•5^5^sf:!f5f:if3<■>^:f:f5f•:f:f5f5f>^•:f>^•5{•i^:^•5f•5f5f•^■5^5t•*3f•<■:f5f^**5^:f5^:f3f■3f5^/ 


VVllLC^  f  ) 

write(’^/)'0.5m  frost  will  prevent  breakthrough  of  most  heavy' 
write(*/)'equipement  but  mobility  may  be  limited  by  crushing' 
write(*/)'failure  of  the  surface  material.' 
write(^ 
end  if 


return 

end 
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APPENDIX  B.  FORTRAN  CODE  USING  NRMM/CAMMS  VARIABLES  AND  FORMAT 


Variables 


SCENARIO  (TERRAIN) 


Figure  Bl.  Cold  regions  variable 
description  for  NRFAM II  ver.  2.5.0. 


Scenario  variable 
ISMODL 


No 


t 
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Figure  B2.  Logic  in  the  cold  regions 
module  ofNRMM  II  ver.  2.5.0. 


Clayey? 


Figure  B2  (cant'd).  Logic  in  the  cold  regions  module  ofNRMM  II  ver.  2.5.0. 


C  CRLSNOW  COLD  REGIONS  MODEL  28  March  94 


c 

CCRLSNOW^ 
C  * 


LOGICAL  FUNCTION  CRLSNOW() 

C 

C  18  Jun  92  Implemented  in  CAMMS/NRMM  by  R.B.Ahlvin 
C  15  Mar  94  Paul  Richmond 
C  28  Mar  94  Karen  Faran 
C 

C  Inputs: 

C  Common  /DERIVE/ 

C  JPSI  Current  deflection  scenario  index 

C  Common  /PREP/ 

C  GCA(NAMBLXNJPSI)  Ground  contact  area  for  each  traction  element 
C  Common /TERRAN/ 
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C  DSNOW  Depth  of  snow  cover  (0=no  snow)  [in] 

C  DTHAW  Depth  of  thawing  [in] 

C  ITSURF  Hard  surface  condition:  [l]=dry,  2=wet,  3=ice-covered 

C  SIGMA  Density  of  snow  [gm/cm''3) 

C  Common /VEHICL/ 

C  CL  Minimum  ground  clearance  [in] 

C  DIAW (NAMELY)  Undeflected  tire  diameter,  each  assembly  [in] 

C  NAMELY  #traction  element  assemblies  (axles  and  tracks) 

C  NVEH(NAMELY)  Vehicle  traction  assembly  type:  0=tracked,  l=wheeled 
C  NWHL(NAMELY)  #tires  on  each  wheeled  assembly 
C  SECTW(NAMELY)  Tire  nominal  section  width  [in] 

C  TRAKLN(NAMELY)  Length  of  track  on  ground  [in] 

C  TRAKWD(NAMELY)  Track  width  (one  side)  for  each  tracked  assembly  [in] 
C  WGHT (NAMELY)  Weight  on  each  vehicle  assembly  [lbs] 

C  WT(NAMELY)  Tread  widths,  each  assembly  (center  to  center)  [in] 

C  WTE(NAMELY)  Min.  width  between  traction  elements  [in] 

C 

C  Outputs: 

C  Common  /DERIVE/ 

C  DOWPE(NAMELY)  Drawbar  coefficient  for  each  assembly 

C  RTOWPB(NAMBLY)  Motion  res.  coef.,  powered  or  braked,  each  element 
C  RTOWT(NAMBLY)  Motion  res.  coef.,  towed,  each  element 
C  Common  /TPREP/ 

C  IBASE  CRREL  model  base  code: 

C  0  =  unassigned 

C  1  =  super-highway,  primary  road 

C  2  =  secondary  road 

C  3  =  off-road  &  trail 

C  4  =  Ice 

C  5  =  soft  soil 

C 

C  External: 

C  RHARD  Function  to  return  hard  surface  resistance  for  snow 
C 

C  Internal: 

C  AA  Track  approach  angle  (rad) 

C  AREAl  Contact  area  single  traction  element  (in'^2) 

C  DEFSEC  1  minus-the  deflected  radius  divided  by  the  section  width  (in) 

C  FACT  Conversion  factor  for  N/cm'^2  to  psi 
C  FLAG  Deep  snow  multiplier 
C  I  Traction  assembly  index 

C  IBASE  Terrain  type  flag  as  follows: 

C  1  =  road 

C  2  =  frozen  or  strong  soil 

C  3  =  unfrozen  soil  (RCI<100) 

C  4  =  ice 

C  PMAX  Maximum  suspension  element  contact  pressure  (psi) 

C  PN  Number  of  vehicle  passes 

C  PRIKPA  Contact  pressure  of  single  traction  element  Kilo-Pascals 

C  PRESl  Contact  pressure  of  single  traction  element  (psi) 

C  RADSI  Radius  of  current  wheeled  traction  element  (in) 

C  RRESl  Motion  resistance  from  snow  (lb) 

C  RTOW  Temporary  resistance  coefficient 
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c 

SIGMAO 

Input  snow  specific  gravity 

c 

SIGMAF 

Final  snow  specific  gravity 

c 

SLIP 

Vehicle  wheel  slip 

c 

SOILZ 

Sinkage  in  soil  of  current  traction  element  (in) 

c 

SUMZ 

Total  vehicle  sinkage 

c 

TFOW 

Temporary  traction  coefficient 

c 

TRCNl 

Maximum  traction  for  i-th  tire  (lb) 

c 

XL2 

Snow  contact  length  for  i-th  suspension  element  (in) 

c 

XNWORT 

Number  of  wheels  or  tracks  on  suspension  assembly 

c 

XWIDl 

Width  of  one  traction  element  (in) 

c 

WGHTl 

Weight  on  single  traction  element  (wheel  or  track)  (lb) 

c 

X2KPA 

Conversion  factor  from  psi  to  KPA 

c 

ZMAX 

Maximum  sinkage  in  snow  (in) 

c 

IMPLICIT  NONE 
INCLUDE  'nrmmdefs.inc' 

INCLUDE  'nrmmdrvc.inc' 

INCLUDE  'nrmmcntc.inc' 

INCLUDE  'nrmmvehc.inc' 

INCLUDE  'nrmmvppc.inc' 

INCLUDE  'nrmmscnc.inc' 

INCLUDE  'nrmmterc.inc' 

INCLUDE  'nrmmtppc.inc' 

C 

INTEGER  I 
REALACOS 

REAL  AA,  AREAl,  DEFSEQ  FLAG 
REAL  PMAX,  PN,  PRESl,  RADSI,  RRESl,  RTOW 
REAL  SIGMAO,  SIGMAF,  SLIP,  SOILZ,  SUMZ 
REAL  TFOW,  TRCNl,  WGHTl,  X2KPA,  XL2 
REAL  XNWORT,  XWIDl,  ZMAX 
C 

REAL  RHARD 
EXTERNAL  RHARD 
C 

C  Conversion  factor  for  psi  to  KPA 

PARAMETER  ( X2KPA  =  XM2IN  *  XM2IN  /  1000.0  /  XN2LB ) 

C 

C  Diagnostic  printout  of  input  data: 

IF(  KIV(3)  .GT.  YES  )THEN 

CALL  DIAGIN('!0’^**Routine  CRLSNOW(IV3)  7/  'CRREL  snow  model 
CALLDIAG('CL  DIAW  DSNOW  DTHAW  GCA'  ) 
CALLDIAG('ITSURFITUT  NAMELY  NVEH  NWHU  ) 

CALL  DIAG(  'SECTW  SIGMA  TRAKLN  TRAKWD  WGHT' ) 

CALL  DIAG(  'WT  WTE' ) 

END  IF 
C 

CRLSNOW  =  .TRUE. 

C 

C  ...  Set  Snow  model  flag,  IBASE 
IF(  ITSURF.EQ.3  )THEN 
C  Ice 

IBASE  =  4 

ELSE  IF(  ITUT.GE.11.AND.ITUT.LE.13  )THEN 


42 


uu  uuuuu  uuu 


C  Hard  surface 
IBASE  =  1 

C  ELSE  IF(  DTHAW  .LE.  0.0  )THEN 

ELSE  IF(  DFREEZ  .GT.  0.0  .AND.  DTHAW  .LE.  0.0  )THEN 
C  Frozen  soil  surface 
IBASE  =  2 
ELSE 

C  LFnfrozen  soil 
IBASE  =  3 
END  IF 
C 

SIGMAO  =  SIGMA 

...  Find  maximum  contact  pressure 
PMAX  =  0.0 
DO  10  1=1,  NAMBLY 
IF(  NVEH(I)  .EQ.  WHEEL  )THEN 
AREAl  =  GCA(LJPSI) 

XNWORT  =  NWHL(I) 

ELSE 

AREAl  =  GCA(LJPSI)  /  2.0 
XNWORT  =  2.0 
END  IF 

Weight  on  single  traction  element 
WGHTl  =  WGHT(I)  /  XNWORT 
Maximum  contact  pressure  of  all  traction  elements 
PMAX  =  AMAX1(  PMAX,  WGHTl /AREAl ) 

10  CONTINUE 

...  Compute  final  density  sigmaf 
IF(  PMAX  .GT.  101.5  )  THEN 
SIGMAF  =  0.65 

ELSE  IF(  PMAX  .GT.  50.8  )  THEN 
SIGMAF  =  0.60 

ELSE  IF(  PMAX  .GT.  30.5  )  THEN 
SIGMAF  =  0.55 
ELSE 

SIGMAF  =  0.5 
END  IF 

FLAG  =  1.0 
SUMZ  =  0.0 

Initial  density  and  final  density  information 
ZMAX  =  AMAX1(  0.0,  DSNOW  *  (1.  -  (SIGMAO /SIGMAF)) ) 

IF(KIV(3).GT.YES)WRITE(LUN1,601)PMAX,SIGMA0,SIGMAF,ZMAX 
601  FORMAT(' Max  elmnt  cntct  pres  PMAX  =',F8.4/ 

&  '  Inp  snow  density  SIGMAO  =',F8.4/ 

&  '  Final  snow  density  SIGMAF  =',F8.4/ 

&  '  Max  avail  sinkage  ZMAX  =',F8.4) 

C 

DO  20  I  =  1,  NAMBLY 
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Compute  terrain  rolling  resistance 

(based  on  first  assembly  characteristics) 

IF(  NVEH(l)  .EQ.  WHEEL  )THEN 
RADSI  =  DlAW(l)/2.0 
XWIDl  -  SECTW(l) 

DEFSEC  =  (1.0-(DFLCT(LNJPS1)/SECTW(1))) 
ELSE 

XWIDl  =  TRAKWD(l) 

END  IF 

Weight  on  a  single  traction  element 
WGHTl  =  WGHT(l)  /  XNWORT 
PRESl  =  WGHTl /AREAl 


Account  for  weak  soil  base 

number  of  passes  and  slip  assumed  constant  for  following  eqns. 
Number  of  vehicle  passes  made 
PN  =  1.0 
Wheel  slip 
SLIP  -  0.05 

IF(  NVEH(I)  .EQ.  WHEEL  )  THEN 

IF  (IBASE  .EQ.  3)  THEN 
Unfrozen  soil 

SOILZ  =  (10.*  RADSrpn**0.5)/((CIRCI*XWIDl*RADSI*2.0)/ 
(WGHTl  *  DEFSEC**(3./2.)  * 

SLIP**(l./5.)))**(5.0/3.0) 

IF  (SOILZ  .LT.  0.0)  SOILZ  =  0.0 
ELSE 

SOILZ  =  0.0 
END  IF 

SUMZ  =  ZMAX  +  SOILZ 

XL2  =  RADSI  *  ACOS((RADSI-SUMZ)/RADSI) 


Account  for  tire  curvature/ deep  snow  flag 
IF(  (SIGMAO.LT.0.15)  .AND.  (SUMZ.GE.RADSI) )  FLAG=1.5 
IF(  SIGMAO  .GE.  0.15  )  THEN 
IF(  SUMZ  .GE.  (RADSF2.0/3.0) )  FLAG  =  1.5 
IF(  SUMZ  .GE.  RADSI )  FLAG  =  2.5 
END  IF 
Belly  Drag 

IF(  SUMZ  .GE.  CL)  FLAG  =  FLAG  +  1.5 


Tracked 

IF  (IBASE  .EQ.  3)  THEN 
SOILZ  =  TRAKLN(l)  *  pn**0.5  *  0.00443  * 

&  EXP((5.889  *  WGHTl)  /  (CIRCI  *  XWIDl  *  TRAKLN(l)))  - 

ELSE 
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SOILZ  =  0.0 
END  IF 

c 

SUMZ  =  ZMAX  +  SOILZ 

IF(KIV(3).GT.YES)WRITE(LUN1,602)SOILZ 
602  FORMAT('  Sinkage  in  soil  SOILZ  =',F8.4) 

Since  NRMM  does  not  have  the  approach  angle  of  a  track  available 
in  the  input  data,  is  assumed  to  be  about  26  deg.  for 
all  tracked  vehicles. 

AA  =  26.0  *  DEGRAD 
XL2  =  (  SUMZ)  /  SIN(  AA ) 

Belly  drag ,  low  (flag=2.0)  or  high  {flag=4.0)  density 

IF(  SUMZ  .GE.  CL  )  THEN 
IF(  SIGMAO  .LT.  0.15)  THEN 
FLAG  =  2.5 
ELSE 

FLAG  =  4.0 
END  IF 
END  IF 
END  IF 

Warning  for  deep  snow 
CALL  ERROR(  'CRLSNO(IV3)', 

&  '***Waming***  Vehicle  belly  dragging  in  deep  snow') 


....  Calculate  terrain  resistance  for  up  to  4  units  for  wheeled  vehicles 

IF  (NVEH(I).EQ.  WHEEL)  THEN 
IF  (LLE.4)  THEN 

RRESl  =  ( (  68.1  *  ( 0.10*  SIGMAO  *  (XL2/XCM2IN)  * 

&  (XWID1/XCM2IN)  )**0.914)  *  FLAG  )  *  XN2LB 

ELSE 

RRESl  =  0.0 
END  IF 
END  IF 

1-unit  limit  for  tracks 

IF  (NVEH(I).EQ.  TRACK)  THEN 
IF  (LLT.2)  THEN 

RRESl  =  ( ( 157.54  *  ( 0.10*  SIGMAO  *  (XL2/XCM2IN)  * 

&  (XWIDl  /XCM2IN)  )**0.833)  *  FLAG  )  *  XN2LB 

ELSE 

RRESl  =  0.0 
END  IF 
END  IF 

Terrain  resistance  on  packed  snow  is  zero 
IF  (SIGMAO  .GE.  0.55)  RRESl  =  0.0 
C  Terrain  resistance  on  ice  is  zero 
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IF  ((IBASE.EQ.4)  .AND.  ( DSNOW  .LE.  (1.0/XCM2IN)))  RRESl  =  0.0 


..  Traction  calculation 

IF(  NVEH(I)  .EQ.  WHEEL  )THEN 
AREAl  =GCA(IJPSI) 

XNWORT  =  NWHL(I) 

XWIDl  =  SECTW(I) 

ELSE 

AREAl  =  GCA(IJPSI)  /  2.0 
XNWORT  =  2.0 
XWIDl  =  TRAKWD(I) 

END  IP 

Weight  on  a  single  traction  element 
WGHTl  =  WGHT(I)  /  XNWORT 
PRESl  =  WGHTl /AREAl 

IP(  IBASE  .NE.  4  )THEN 

Snow  on  road  (IBASE=1) 

Snow  on  frozen  or  strong  (RCI>100)  soil  (IBASE=2) 
Snow  on  unfrozen  soil  (RCI<=100)  (IBASE=3) 

IF  (SIGMAO.GE.0.55)  THEN 
hard  packed  snow 

TRCNl  =  0.321  (PRES1*X2KPA)’^*0.97 
ELSE 

TRCNl  =  0.851  *  (PRES1*X2KPA)’^*0.823 
END  IF 
ELSE 

Snow  on  ice  base  (IBASE=4) 

IF  (  DSNOW  .LE.  (1.0/XCM2IN))  THEN 
TRCNl  =  0.1  (PRESl ’^X2KPA) 

ELSE 

IF  (SIGMAO.GE.0.55)  THEN 
Disturbed  snow 

TRCNl  =  0.321  *  (PRES1*X2KPA)*’'0.97 
ELSE 

Undisturbed  snow 

TRCNl  =  0.127  *  (PRESl’^X2KPA)*n.06 
END  IF 
END  IF 
END  IF 


TRCNl  is  a  stress  at  this  point 
TRCNl  =  (TRCNl /X2KPA)  *  AREAl 
C 

TFOW(I)  =  TRCNl  *  XNWORT  /  WGHT(I) 

RTOW  =  RRESl  *  XNWORT  /  WGHT(I)  +  RHARD(I) 
C 

RTOWPB(I)=  RTOW 
RTOWT(I)  =  RTOW 
DOWPB(I)  =  TFOW(I)  -  RTOW 
C 

20  CONTINUE 
C 
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99  RETURN 
END 

C  RHARD  Hard  surface  resistance  for  snow  19  Jun  92 

c 

CRHARD* 

C 

REAL  FUNCTION  RHARD(I) 

C 

C  This  routine  computes  the  motion  resistance  for  the  hard  surfaces 
C  to  be  added  to  the  resistance  computed  for  snow. 

C 

C  19  Jun  92  Origional  edit  (taken  from  part  of  origional  NRMM  snow  model) 
C 


C  Inputs: 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


I 

Common  /PREP/ 
CPFFG(NAMBLY) 
Common  /TPREP/ 
ITUT 


Suspension  assembly  index 
F-G  contact  pressure  factor  for  each  element 


Terrain  type  code: 

1  =  off-road  (areal)terrain 

2  =  water-covered  terrain 

11  =  on-road  super-highway 

12  =  on-road  primary  road 

13  =  on-road  secondary  road 

14  =  on-road  trail 
Common  /VEHICL/ 

NVEH(NAMBLY)  Vehicle  traction  assembly  type:  0=tracked,  l=wheeled 


C  Outputs: 

C  RHARD  Hard  surface  resistance  coefficient 
C 

IMPLICIT  NONE 
INCLUDE  'nrmmdefs.inc' 

INCLUDE  'nrmmcntc.inc' 

INCLUDE  'nrmmvehc.inc' 

INCLUDE  'nrmmvppc.inc' 

INCLUDE  'nrmmtppc.inc' 

C 


INTEGER  I 


C 

IF(  NVEH(I)  .EQ.  WHEEL  )THEN 
C  Hard  surface  resistance  coefficients  (wheeled) 
IF(  ITUT.EQ.l  .OR.  ITUT.EQ.14  )THEN 
C  Off-road  and  trails  (no  ice) 

IF(  CPFFG(I).GE.4.0  )THEN 
RHARD  =  0.0175 
ELSE 

RHARD  =  0.015 
END  IF 

ELSE  IF  (ITUT.EQ.il  .OR.  ITUT.EQ.12  )THEN 
C  Super-highways  &  primry  roads  or  ice 
RHARD  =  0.015 


47 


ELSE 

C  Secondary  roads  &  all  others 

RHARD  =  0.025 
END  IF 
ELSE 

C  Hard  surface  resistance  coefficients  (tracked) 

IF(  ITUT.EQ.l  .OR.  ITUT.EQ.14  )THEN 
C  Off-road  and  trails 

RHARD  =  0.0525 

ELSE  IF  (ITUT.EQ.il  .OR.  ITUT.EQ.12  )THEN 
C  Super-highways  &  prhnry  roads 

RHARD  =  0.0375 
ELSE 

C  Secondary  roads  &  all  others 

RHARD  =  0.045 
END  IF 
END  IF 
C 

IF(KIV(3).GT.YES)WRITE(LUN1,601)RHARD 
601  FORMAT('  Hard  Sfc  Resist.  Coef.  RHARD  =',F8.4) 
C 

RETURN 

END 


C  THAW  CRREL  Thawing  ground  traction  model  6  Jan  94 


C 

CTHAW* 

C 


LOGICAL  FUNCTION  THAW() 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


6  May  94  Implemented  in  NRMM II 
Inputs: 

Common  /DERIVE/ 

DOWPB(NAMBLY)  Drawbar  coefficient  for  each  assembly 
LAYER  Soil  layer  index  (Critical  layer):  l=:0"-6",  2=6"-12" 

RTOWPB(NAMBLY)  Motion  res.  coef.,  powered  or  braked,  each  element 
RTOWT(NAMBLY)  Motion  res.  coef.,  towed,  each  element 
Common  /TERRAN/ 

DTHAW  Ground  thawing  depth  [in] 

DSNOW  Snow  depTH  [in] 

Common  /VEHICL/ 

NAMBLY  #traction  element  assemblies  (axles  and  tracks) 


Outputs: 

THAW  Return  status  .TRUE.  =  O-K 

Coirunon  /DERIVE/ 

DOWPB(NAMBLY) 

MRMUL(NAMBLY) 

RTOWPB(NAMBLY) 

RTOWT(NAMBLY) 

TFMUL(NAMBLY) 

TFOW(NAMBLY) 


Drawbar  coefficient  for  each  assembly 
Winter  model  soil  resistance  multiplier,  each  element 
Motion  res.  coef.,  powered  or  braked,  each  element 
Motion  res.  coef.,  towed,  each  element 
Winter  model  traction  multiplier,  each  element 
Traction  coefficient  from  soil,  each  element 
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c 

C  External: 

C  MRMULT  Soil  motion  coefficient  multiplier  for  thawing  soils 
C  TRMULT  Soil  traction  coefficient  multiplier  for  thawing  soils 
C 

C  Internal; 

C  I  Suspension  assembly  index 
C  STATOK  Temporary  status  flag 
C 

IMPLICIT  NONE 
INCLUDE  'nrmmdefs.inc' 

INCLUDE  'nrmmdrvc.inc' 

INCLUDE  'nrmmcntc.inc' 

INCLUDE  'nrmmvehc.inc' 

INCLUDE  'nrmmvppc.mc' 

INCLUDE  'nrmmterc.inc' 

INCLUDE  'nrmmtppc.inc' 

C 

LOGICAL  STATOK 
INTEGER  I 

REAL  MRMULT,  TRMULT 
EXTERNAL  MRMULT,  TRMULT 


Diagnostic  printout  of  input  data: 

IF(  KIV(3)  .GT.  YES  )THEN 

CALL  DIAGIN('!0***Routine  THAW(IV3) '// 

&  'CRREL  Thawing  soil  traction  model') 

CALL  DIAG(  'DOWPB  DTHAW  DSNOW  LAYER  NAMELY' ) 
CALL  DIAG(  'RTOWPB  RTOWT' ) 

END  IF 

STATOK=.TRUE. 

DO  10 1=1,  NAMELY 

TFMUL(I)  =  TRMULTO 
IF(STATOK)STATOK  =  TFMUL(I).GT.O 
TFOW(I)  =  (DOWPB(I)+RTOWPB(I))*TFMUL(I) 

MRMUL(I)  =  MRMULT(I) 

IF(STATOK)STATOK  =  MRMUL(I).LT.999. 

RTOWPB(I)  =  RTOWPB(I)*MRMUL{I) 

RTOWT(I)  =  RTOWT(I)*MRMUL(I) 

DOWPB(I)  =  TFOW(I)  -  RTOWPB(I) 

10  CONTINUE 
THAW  =  STATOK 


Check  for  validity 
IF(  LAYER  .EQ.  2  ) 

&  CALL  ERROR(  'THAW','Critical  depth  >  6"' ) 
IF(  DSNOW  .GT.  1./XCM2IN) 

&  CALL  ERROR(  'THAW','Snow  cover  >  1  cm.' ) 


RETURN 

END 


C  TRMULT  Traction  multiplier  for  thawing  soils  6  Jan  94 
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c 

CTRMULT=^ 

C 

^  :f  St  >(■  5f  >e- 5t  5(- 

REAL  FUNCTION  TRMULT() 

C 

C  May  92  Model  based  on  US  ARMY  CRREL  research  done  by  Sally  Shoop 
C  Coded  by  Jim  Slota,  OptiMetrics  for  CRREL.  (SSM  6.2) 

C  23  Jun  92  Implemented  in  CAMMS  by  R.  B.  Ahlvin,  WES 
C  6  Jan  94  Implemented  in  NRMM II  by  R.B.Ahlvin 
C 

C  This  function  returns  a  multiplier  value  to  modify  drawbar-pull 
C  in  lbs.  Traction  is  reduced  for  travel  over  a  thawing  soil  layer 
C  that  overlays  a  frozen  soil  layer.  Calculations  are  applicable 
C  only  for  sandy  soils  with  little  or  no  snow  cover,  and  for 
C  vehicles  with  a  critical  depth  of  less  than  6  inches.  The  wetness 
C  index,  as  calculated  by  the  SMSP  model,  is  also  required. 

C 


C  Inputs: 

C  Common  /TERRAN/ 

DTHAW  Depth  of  thawing  is  soil  [in] 

KUSCS  uses  soil  type  code:  1=SW,  2=SP,  3=SM,  4=SC,  5=SMSC, 

6=CL,  7=ML,  8=CLML,  9=CH,  10=MH,  ll=OL,  12=OH,  13=WATER, 
14=PAVEMENT,  15=ROCK,  16=GW,  17-GP,  18-GM,  19=GC,  20=Pt 
SMSP  wetness  index:  0=arid,  l=dry,  2=average,  3=wet, 
4=saturated,  5=waterlogged 


C 
C 

c 
c 
c 
c 

C  Output: 


KWI 


C  TRMULT  Tractio  multiplier 
C 


IMPLICIT  NONE 
INCLUDE  'nrmmdefs.inc' 
INCLUDE  'nrmmcntc.inc' 
INCLUDE  'nrmmterc.inc' 
INCLUDE  'nrmmvehc.inc' 


C 


INTEGER  I,J,  LENSTR 
REAL  MRMULT,  THAWCM 
LOGICAL  SVALID(20) 

CHARACTER  *4  ZUSCS(20) 

CHARACTER  *5  ZVTYPE 
EXTERNAL  LENSTR 

DATAZUSCS/'SW','SP','SM  ','SC  ','SMSC', 
&  'CL','ML','CLMU,'CH  ','MH  ', 

&  'OU,'OH','WATR',TVMT','ROCK', 

&  'GW','GP','GM  ','GC  ',Tt  7 

C 

DATA  (SVALID(J),J=  1,  5) 

&  /.TRUE.  ,.TRUE.  ,.TRUE.,  .TRUE.,  .TRUE./ 
C  SW  SP  SM  SC  SMSC 

DATA  (SVALID(J),J=  6,10) 

&  /.TRUE.  ,.FALSE.,.FALSE.,.FALSE.,.FALSE./ 
C  CL  ML  CLML  CH  MH 

DATA  (SVALID(J),J=11,15) 

&  /.FALSE.,.FALSE.,.FALSE.,.FALSE.,.FALSE./ 


50 


C  OL  OH  WATR  PVMT  ROCK 

DATA  (SVALIDQ)J=16,20) 

&  /.TRUE.  ,.TRUE.  ,.TRUE.  ,.TRUE.  ,.TRUE./ 

C  GW  GP  GM  GC  Pt 

C 

THAWCM  =  DTHAW  /  XCM2IN 

IF(  (THAWCM  .GT.  2.5)  .AND.  (THAWCM  .LT.  15.0))  THEN 
TRMULT  =  2.379  *  (1.0/(THAWCM*THAWCM))  +  0.619 
ELSE  IF  (THAWCM  .GE.  15.0)  THEN 
TRMULT  =  0.63 
ELSE 

TRMULT  =  1.0 
END  IF 

IF(  .NOT.SVALID(KUSCS) )  THEN 

WRITE(ERRMSG,9901)KUSCS,ZUSCS(KUSCS)(1:LENSTR(ZUSCS(KUSCS))) 
9901  FORMAT('Relations  not  implemented  for 
&  'soil  type  KUSCS=',I2/  (',A/)': 

&  '  NVEH(I)='T1/  (',A/)') 

CALL  ERROR(  'TRMULT',  ERRMSG  ) 

END  IF 


RETURN 


C  MRMULT  Motion  resistance  multiplier  for  thawing  soils  23  Jun  92 


C 

CMRMULT* 

C 

^  X- >(■  X- ♦  if- St- if- X- if 

ENTRY  MRMULT(I) 

C 

C  May  92  Model  based  on  US  ARMY  CRREL  research  done  by  Sally  Shoop 
C  Coded  by  Jim  Slota,  OptiMetrics  for  CRREL.  (SSM  6.2) 

C  6  May  94  Implemented  in  NRMM  II  by  R.B.AhIvin 
C 

C  This  function  returns  a  multiplier  value  to  motion  resistance 
C  values.  This  applies  to  travel  over  a  thawing  soil  layer 
C  overlaying  a  frozen  layer.  Calculations  are  applicable  only  for 
C  wheeled  vehicles  on  sandy  soils  with  little  or  no  sno  cover, 

C  and  for  vehicles  with  a  critical  depth  of  less  than  6  inches.  The 
C  wetness  index,  as  calculated  by  the  SMSP  model,  is  also  required. 

C 


C  Inputs: 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


I  Suspension  assembly  index 

Common  /TERRAN/ 

DTHAW  Depth  of  thawing  is  soil  [in] 

KUSCS  uses  soil  type  code:  1=SW,  2=SP,  3=SM,  4=SC,  5=SMSC, 

6=CL,  7=ML,  8=CLML,  9=CH,  10=MH,  ll=OL,  12=OH,  13= WATER, 
14=PAVEMENT,  15=ROCK,  16=GW,  17=GP,  18=GM,  19=GC,  20=Pt 
KWI  SMSP  wetness  index:  0=arid,  l=dry,  2=average,  3=wet, 

4=saturated,  5=waterlogged 
Common  /VEHICL/ 

NVEH(NAMBLY)  Vehicle  traction  assembly  type:  0=tracked,  l=wheeled 
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C  Output: 

C  MRMULT  Motion  resistance  multiplier 
C 

MRMULT  =1.0 
C 

THAWCM  =  DTHAW  /  XCM2IN 

C  The  following  3  lines  handle  cases  for  which  THAWCM  <  2.5  cm 
IF(  THAWCM  .LE.  2.5  )THEN 
MRMULT  =  1.0 

ELSE  IF(  (THAWCM  .GT.  2.5)  .AND.  (THAWCM  .LT.  4.0))  THEN 
MRMULT  =  2.883  *  THAWCM  -  6.056 

ELSE  IF(  (THAWCM  .GE.  4.0)  .AND.  (THAWCM  .LT.  8.0))  THEN 
MRMULT  =  -0.22*THAWCM*THAWCM  +  3.54*THAWCM  -  5.24 
ELSE  IF(  (THAWCM  .GE.  8.)  .AND.  (THAWCM  .LT.  12.0))  THEN 
MRMULT  =  0.225*THAWCM  +  7.2167 
ELSE 

MRMULT  =  10.0 
END  IF 

IF(  NVEH(I).EQ.  TRACK  ) 

&  CALL  ERROR( 'MRMULT', 

&  'Multiplier  may  not  be  valid  for  TRACKED  elements' ) 

IE(  .NOT.SVALID(KUSCS) )  THEN 

WRITE(ERRMSG,9901)KUSCS,ZUSCS(KUSCS)(1:LENSTR(ZUSCS(KUSCS))), 
&  NVEH(I),ZVTYPE 

CALL  ERROR(  'MRMULT',  ERRMSG  ) 

END  IF 
C 

RETURN 

END 
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APPENDIX  C.  TRACTION  COEFFICIENTS  ON  PACKED  SNOW 


Traction  on  groomed,  shallow  snow  was  measured  with  three  different  vehicles —  the  CRREL  Instru¬ 
mented  Vehicle  (CIV),  the  Saab  Friction  Tester  and  the  Unrroyal-Goodrich  Traction  Tester  (U-G) — in  Shoop 
(1993a).  Several  methods  of  reporting  the  traction  coefficient  (p)  were  used: 

•  The  maximum  value  (mnax)i 

•  An  average  peak  value  (|ipeak)i 

•  The  value  at  12%  slip  (Pi2%slip)i 

•  The  SAE-specified  p  (Psae)- 

These  measurements  are  compared  to  other  published  traction  coefficients  on  groomed  snow  in  Figure 
Cl.  The  predicted  traction  values,  using  CRREL's  Cold  Regions  Mobility  Models,  is  expressed  in  coeffi¬ 
cient  form  and  is  shown  by  the  open  circles  at  the  bottom  of  the  figure. 


Snow 


•  pciv  Max 
•  PciV  Ave.  Peak 
Pciv  12%  Slip 
Pciv  sae 


hsaab12%Slip 


PcivAve.  Peak  (gross) 


•  hu-G 

Psaab12%Slip 

•  Pu-G  SAE 

•  Pciv  Max 

•  Pu-oMax 
•  hcivAve.  Peak 

•  Pciv  12%  Slip 

•  hu-G^AE 
I  PcivSAE  ; 

Pciv  Ave.  Peak  (gross) 

Pciv  Ave.  Peak  (net) 

■  Pciv  Ave.  Peak  (gross) 
p  Packed  Snow  (SME,  1973) 


•  Data  from  Shoop  1993 
■  Data  from  Richmond  et  al.  1990,  SME  1977 
O  Predicted  Gross  Traction 


Snow  Road 
>■  2.5  cm,  0.25  g/cm® 
Snow  over  Hard  Pack 


Snow  Test  Section  '92 
10  cm,  0.54  g/cm® 
Snow  over  Hard  Pack 
(Tair  =  -2' 


P  f  _ -iC 

'  snow  “  ' 


C) 


Snow  Test  Section  '93 
2  cm,  0.59  g/cm® 
Snow  over  Hard  Pack 

(Tair  =  -0.6° 


C.  ''‘snQW  =  ~2.8° 


C) 


6  cm, 0.59  g/cm^  Groomed 
and  Packed  Snow 
(Richmond  etai.,  1990). 


O  Psaab  Predicted.  Hard  Packed  Snow 
O  Pciv  Predicted,  Hard  Packed  Snow 
O  Psaab  Predicted,  Undisturbed  Snow 

O  Pciv  Predicted,  Undisturbed  Snow 


0.2 


0.4 


0.6 


Traction  Coefficient  (p) 

Figure  Cl.  Comparison  of  measured,  published  and  predicted  values  of  traction  on  shallow 
groomed  snow. 
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APPENDIX  D:  NRMM  CHECKOUT  DATA 


TERRAIN  UNIT  INPUT 

38  8  Checkout  terrain  for  CRREL  winter  effects  model  6  May  94 

! 

!  1)  Best  possible  on-road  terrain 

!  NTU  TUID  AREA  IROAD ISCOND  GRADE  ACTRMS  RDA  EANG  RADC 
l'BstON'1.00  1  1  0  0.0  999  0  5000 

! 

!  2)  Best  possible  off-road  terrain 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

2  'BstCC'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

1 

!  3)  Soft  soil  effects.  Soil  strength  =  100 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

3  'CIlOO'  1.00  0  1  0  0.0  999  SM  2  100  100  99 

I 

!  4)  Soft  soil  effects,  Soil  strength  =  50 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
4'CI50'  1.00  0  1  0  0.0  999  SM  3  50  50  99 

I 

!  5)  Soft  soil  effects.  Soil  strength  =  20 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

5  'CI20'  1.00  0  1  0  0.0  999  SM  4  20  20  99 

! 

!  6)  Frozen  ground,  w/Ice  base 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

6  Tee'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  99.9  0.0  0.0  1 
! 

!  7)  1-in  snow  on  frozen  ground,  w/Ice  base 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
7 'lin/I' 1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  1.0  0.2  99.9  0.0  0.0  1 

I 

!  8)  10-in  snow  on  frozen  ground,  w/Ice  base 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

8  'lOinF  1.00  0  1  0  0.0  999  SM  1  300  300  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  10.0  0.2  99.9  0.0  0.0  1 

I 

!  9)  1-in  snow  on  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

9  'lin'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  1.0  0.2  99.9  0.0  0.0  0 


!  10)  10-in  snow  on  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
10  'lOin'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  99.9  0.0  0.0  0 


!  11)  20-in  snow  on  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
11  '20in'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  20.0  0.2  99.9  0.0  0.0  0 


!  12)  10-in  snow  on  frozen  ground,  density=  0.1 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
12  '10".l'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.1  99.9  0.0  0.0  0 
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!  13)  10-in  snow  on  frozen  ground,  density=  0.3 

!  NTU  TUID  AREA  IROAD ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
13 '10".3' 1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.3  99.9  0.0  0.0  0 


!  14)  10-in  snow  on  frozen  ground,  density=  0.4 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
14  '10”.4r  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.4  99.9  0.0  0.0  0 
! 

!  15)  10-in  snow  on  road 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  BANG  RADC 
15 '10"ON' 1.00  1  5  0  0.0  999  0  5000 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  0.0  0.0  0.0  0 
! 

!  16)  10-in  snow  on  unfrozen  ground,  CI=150 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
16  '10CI150'  1.00  0  1  0  0.0  999  SM  1  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  0.0  0.0  0.0  0 
! 

!  17)  10-in  snow  on  unfrozen  ground,  CI=50 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
17 '10CI50' 1.00  0  1  0  0.0  999  SM  2  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  0.0  0.0  0.0  0 
! 

!  18)  10-in  snow  on  unfrozen  ground,  CI=20 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
18 '10CI20' 1.00  0  1  0  0.0  999  SM  3  20  20  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  0.0  0.0  0.0  0 

I 

!  19)  10-in  snow  on  2-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
19 '2frez' 1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  2.0  0.0  0.0  0 
! 

!  20)  10-in  snow  on  4-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

20  '4frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  10.0  0.2  4.0  0.0  0.0  0 

I 

!  21)  10-in  snow  on  6-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

21  i'6frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  10.0  0.2  6.0  0.0  0.0  0 
! 

!  22)  10-in  snow  on  12-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

22  '12frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  10.0  0.2  12.0  0.0  0.0  0 

I 

!  23)  10-in  snow  on  2-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

23  '2thaw'  1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  10.0  0.2  12.0  2.0  0.0  0 
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!  24)  10-in  snow  on  4-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
24 '4 thaw' 1.00  0  1  0  0.0  999  SM  4  20  20  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  12.0  4.0  0.0  0 
! 

!  25)  10-in  snow  on  8-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
25  '8thaw'  1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  12.0  8.0  0.0  0 

!!  26)  10-in  snow  on  8-in  thawing  ground,  10%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
26 'SthMlO' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  12.0  8.0  10.0  0 


!  27)  10-in  snow  on  8-in  thawing  ground,  30%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
27 '8thM30' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  12.0  8.0  30.0  0 


!  28)  10-in  snow  on  8-in  thawing  ground,  50%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
28 '8thM30' 1.00  0  1  0  0.0  999  SM  4  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  10.0  0.2  12.0  8.0  50.0  0 
! 

!  29)  0-in  snow  on  2-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

29  '2frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  2.0  0.0  0.0  0 
] 

!  30)  0-in  snow  on  4-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

30  '4frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  4.0  0.0  0.0  0 
! 

!  31)  0-in  snow  on  6-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

31  i'6frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  6.0  0.0  0.0  0 
! 

!  32)  0-in  snow  on  12-in  frozen  ground 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

32  '12frez'  1.00  0  1  0  0.0  999  SM  1  300  300  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  0.0  0.0  0 


!  33)  0-in  snow  on  2-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
33 '2thaw' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 

!  34)  0-in  snow  on  4-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
34 '4thaw' 1.00  0  1  0  0.0  999  SM  4  20  20  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  4.0  0.0  0 
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!  35)  0-in  snow  on  8-in  thawing  ground  (using  WI) 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
35 '8thaw' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  8.0  0.0  0 
! 

!  36)  0-in  snow  on  8-in  thawing  ground,  5%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
36 '8thM10' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  8.0  5.0  0 
! 

!  37)  0-in  snow  on  8-in  thawing  ground,  10%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
37 '8thM30' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  8.0  10.0  0 
! 

!  38)  0-in  snow  on  8-in  thawing  ground,  30%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
38  '8thM30'  1.00  0  1  0  0.0  999  SM  4  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  8.0  30.0  0 

54  8  Checkout  terrain  for  CRREL  winter  effects  model  6  May  94 
! 

!  39)  ZTEST#1:2"  frozen  ground,RCIl,2=5,5;KWI=l 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
39'zt-l'l.OO  0  1  0  0.0  999  SM  1  5  5  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  2.0  0.0  0.0  0 
! 

!  40)  ZTEST  #2:2"  frozen  ground,RCIl,2=20,20;KWI=l 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
40'zt-2'1.00  0  1  0  0.0  999  SM  1  20  20  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  2.0  0.0  0.0  0 
1 

!  41)  ZTEST  #3:2"  frozen  ground,RCIl,2=100,100;KWI=l 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
41  'zt-3'  1.00  0  1  0  0.0  999  SM  1  100  100  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  2.0  0.0  0.0  0 

I 

!  42)  ZTEST  #4:2"  frozen  ground,RCIl,2=5,5;  KWI=5 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
42'zt-4'1.00  0  1  0  0.0  999  SM  5  5  5  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  2.0  0.0  0.0  0 
! 

!  43)  ZTEST  #5:2"  frozen  ground,RCIl,2=20,20;KWI=5 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
43'zt-5'1.00  0  1  0  0.0  999  SM  5  20  20  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  2.0  0.0  0.0  0 
! 

!  44)  ZTEST  #6:2"  frozen  ground,RCIl,2=100,100;KWI=5 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

44  'zt-6'  1.00  0  1  0  0.0  999  SM  5  100  100  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  2.0  0.0  0.0  0 
! 

!  45)  0-in  snow  on  2-in  thawing  ground  (using  WI),  wet  conditions 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

45  '2thaw'  1.00  0  1  0  0.0  999  SM  5  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
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!  46)  SOIL  TEST#1:  KWI=5,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  SW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
46'st-l'l.OO  0  1  0  0.0  999  SW  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
1 

!  47)  SOIL  TEST#2:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

47  'st-2'  1.00  0  1  0  0.0  999  SP  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 

I 

!  48)  SOIL  TEST#3:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SM 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

48  'st-3'  1.00  0  1  0  0.0  999  SW  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  49)  SOIL  TEST#4:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
49'st-4'1.00  0  1  0  0.0  999  SC  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
] 

!  50)  SOIL  TEST#5:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SMSC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
50'st-5'1.00  0  1  0  0.0  999  SMSC  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  51)  SOIL  TEST#6:  KWI=5,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  CL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

51  'st-6'  1.00  0  1  0  0.0  999  CL  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 

!  52)  SOIL  TEST#7:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  ML 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

52  'st-7'  1.00  0  1  0  0.0  999  ML  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  53)  SOIL  TEST#8:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CLML 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

53  'st-8'  1.00  0  1  0  0.0  999  CLML  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 

I 

!  54)  SOIL  TEST#9:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

54  'st-9'  1.00  0  1  0  0.0  999  CH  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 

I 

!  55)  SOIL  TEST#10:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  MH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

55  'st-10'  1.00  0  1  0  0.0  999  MH  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 

!  56)  SOIL  TEST#11:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  OL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

56  'st-3'  1.00  0  1  0  0.0  999  OL  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
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!  57)  SOIL  TEST#12;  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  OH 
!  NTU  TUID  AREA  IROAD ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
57  'st-12'  1.00  0  1  0  0.0  999  OH  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  58)  SOIL  TEST#13:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  WATR 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
58 'st-13' 1.00  0  1  0  0.0  999  WATR  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 

I 

!  59)  SOIL  TEST#14:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  PAVE 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
59 'st-14' 1.00  0  1  0  0.0  999  PAVE  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  60)  SOIL  TEST#15:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  ROCK 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
60  'st-15'  1.00  0  1  0  0.0  999  ROCK  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 

!!  61)  SOIL  TEST#16:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
61 'st-16' 1.00  0  1  0  0.0  999  GW  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 

!  62)  SOIL  TEST#17:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

62  'st-17'  1.00  0  1  0  0.0  999  GP  5  150  150  99 

1  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  63)  SOIL  TEST#18:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GM 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

63  'st-18'  1.00  0  1  0  0.0  999  GM  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  64)  SOIL  TEST#19:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
64 'st-19' 1.00  0  1  0  0.0  999  GC  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  65)  SOIL  TEST#20:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  Pt 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

65  'st-20'  1.00  0  1  0  0.0  999  Pt  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  0.0  0 

!  66)  #2SOIL  TEST#1:  KWHO,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

66  '2st-l'  1.00  0  1  0  0.0  999  SW  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 

!  67)  #2SOIL  TEST#2:  KWHO,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

67  '2st-2'  1.00  0  1  0  0.0  999  SP  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
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!  68)  #2SOIL  TEST#3:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SM 
!  NTU  TUID  AREA  IROAD ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

68  '2st-3'  1.00  0  1  0  0.0  999  SW  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  69)  #2SOIL  TEST#4:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

69  '2st-4'  1.00  0  1  0  0.0  999  SC  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
1 

!  70)  #2SOIL  TEST#5:  KWI=0,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  SMSC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

70  '2st-5'  1.00  0  1  0  0.0  999  SMSC  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  71)  #2SOIL  TEST#6:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

71  '2st-6'  1.00  0  1  0  0.0  999  CL  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
[ 

!  72)  #2SOIL  TEST#7;  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  ML 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
72 '2st-7' 1.00  0  1  0  0.0  999  ML  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

!!  73)  #2SOIL  TEST#8:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CLML 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
73  '2st-8'  1.00  0  1  0  0.0  999  CLML  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  74)  #2SOIL  TEST#9:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
74 '2st-9' 1.00  0  1  0  0.0  999  CH  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  75)  #2SOIL  TEST#10:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  MH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

75  '2st-10'  1.00  0  1  0  0.0  999  MH  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  76)  #2SOIL  TEST#11:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  OL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

76  '2st-3'  1.00  0  1  0  0.0  999  OL  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  77)  #2SOIL  TEST#12:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  OH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

77  '2st-12'  1.00  0  1  0  0.0  999  OH  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  78)  #2SOIL  TEST#13:  KW1=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  WATR 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

78  '2st-13'  1.00  0  1  0  0.0  999  WATR  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  12.0  2.0  40.0  0 
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!  79)  #2SOIL  TEST#14:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  PAVE 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
79  '2st-14'  1.00  0  1  0  0.0  999  PAVE  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  80)  #2SOIL  TEST#15:  KWI=0,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  ROCK 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
80 '2st-15' 1.00  0  1  0  0.0  999  ROCK  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 
[ 

!  81)  #2SOIL  TEST#16:  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
81  '2st-16'  1.00  0  1  0  0.0  999  GW  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  82)  #2SOIL  TEST#17;  KWI=0,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
82 '2st-17' 1.00  0  1  0  0.0  999  GP  0  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  83)  #2SOIL  TEST#18:  KWI=0,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  GM 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
83  '2st-18'  1.00  0  1  0  0.0  999  GM  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

I 

!  84)  #2SOIL  TEST#19:  KWI=0,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  GC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
84 '2st-19' 1.00  0  1  0  0.0  999  GC  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 

!!  85)  #2SOIL  TEST#20:  KWI=0,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  Pt 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
85  '2st-20'  1.00  0  1  0  0.0  999  Pt  0  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  40.0  0 
! 

!  86)  0-in  snow  on  4-in  thawing  ground,  5%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
86 '4thM5' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  4.0  5.0  0 
1 

!  87)  0-in  snow  on  4-in  thawing  ground,  10%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
87 '4thM10' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  4.0  10.0  0 
! 

!  88)  0-in  snow  on  4-in  thawing  ground,  30%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
88  '4thM30'  1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  4.0  30.0  0 
1 

!  89)  0-in  snow  on  4-in  thawing  ground,  50%  moisture  content 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
89 '4thM10' 1.00  0  1  0  0.0  999  SM  3  50  50  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  4.0  50.0  0 
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!  90)  0-in  snow  on  2-in  thawing  ground,  RCI1,2=5 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
90 '2thR5' 1.00  0  1  0  0.0  999  SM  3  5  5  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  91)  0-in  snow  on  2-in  thawing  ground,  RCI1,2=20 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
91  '2thR20'  1.00  0  1  0  0.0  999  SM  3  20  20  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
! 

!  92)  0-in  snow  on  2-in  thawing  ground,  RCI1,2=150 

!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
92 '2thR150' 1.00  0  1  0  0.0  999  SM  3  150  10  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  12.0  2.0  0.0  0 
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! 

!  93)  SOIL  TEST#1:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
93'st-l'l.OO  0  1  0  0.0  999  SW  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!  94)  SOIL  TEST#2:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

94  'st-2'  1.00  0  1  0  0.0  999  SP  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!  95)  SOIL  TEST#3:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SM 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

95  'st-3'  1.00  0  1  0  0.0  999  SW  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

!  96)  SOIL  TEST#4:  KWHS,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  SC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

96  'st-4'  1.00  0  1  0  0.0  999  SC  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 
] 

!  97)  SOIL  TEST#5;  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  SMSC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

97  'st-5'  1.00  0  1  0  0.0  999  SMSC  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!  98)  SOIL  TEST#6:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
98'st-6'1.00  0  1  0  0.0  999  CL  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!  99)  SOIL  TEST#7:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  ML 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
99  'st-7'  1.00  0  1  0  0.0  999  ML  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
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!100)  SOIL  TEST#8:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CLML 
!  NTU  TUID  AREA  IROAD ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

100  'st-8'  1.00  0  1  0  0.0  999  CLML  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!101)  SOIL  TEST#9:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  CH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

101  'st-9'  1.00  0  1  0  0.0  999  CH  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!102)  SOIL  TEST#10:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  MH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
102 'st-10' 1.00  0  1  0  0.0  999  MH  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!103)  SOIL  TEST#11:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  OL 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
103 'st-3' 1.00  0  1  0  0.0  999  OL  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!104)  SOIL  TEST#12:  KWI=5,  RCIl, 2=150,  DFREEZE=12in;DTHAW=2in;  OH 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

104  'st-12'  1.00  0  1  0  0.0  999  OH  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!105)  SOIL  TEST#13:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  WATR 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

105  'st-13'  1.00  0  1  0  0.0  999  WATR  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!106)  SOIL  TEST#14:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  PAVE 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

106  'st-14'  1.00  0  1  0  0.0  999  PAVE  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!107)  SOIL  TEST#15:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  ROCK 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 

107  'st-15'  1.00  0  1  0  0.0  999  ROCK  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 

WIN  0.0  0.2  0.0  0.0  0.0  0 

!108)  SOIL  TEST#16:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GW 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
108 'st-16' 1.00  0  1  0  0.0  999  GW  5  150  150  99 

!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
! 

!109)  SOIL  TEST#17:  KWH5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GP 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
109  'st-17'  1.00  0  1  0  0.0  999  GP  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!110)  SOIL  TEST#18:  KWHS,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GM 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
110 'st-18' 1.00  0  1  0  0.0  999  GM  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
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nil)  SOIL  TEST#19:  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  GC 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
111 'st-19' 1.00  0  1  0  0.0  999  GC  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 

I 

!112)  SOIL  TEST#20;  KWI=5,  RCI1,2=150,  DFREEZE=12in;DTHAW=2in;  Ft 
!  NTU  TUID  AREA  IROAD  ISCOND  GRADE  ACTRMS  RDA  USCS  KWI  RCIl  RCI2  DBROCK 
112  'st-20'  1.00  0  1  0  0.0  999  Ft  5  150  150  99 
!  DSNOW  SIGMA  DFREEZ  DTHAW  TMOIST  ICE 
WIN  0.0  0.2  0.0  0.0  0.0  0 
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TERRAIN  UNIT  OUTPUT 

29-JUL-94  11:17:54  NRMMII  version  2.5.0  last  changed  12  May,  1994 
TERRAIN  1  FIIjE=crrelcko.  ter 

38  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  6  Ma 
SCENARIO  1  FILE=scenario.dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  1  FILE=vehicles/mll3 .crl 

ID=  M113  FOR  CRREL  MODEL  CHECKOUT  3  MAY  94 
Snow  model : CRREL 

Weight  [Lb] :  2®  11700  =  23401 

NTU/  - Gross- (soil  only)  -  [Lb] - - Total- [Lb]  -  Speed 


IBASE  Traction  Resist,  drawbar  DB-Coef. 


Snow  depth=  10 . 0 

in  Density=  0.1 

g/cm''3 

1  0 

9360.1 

9360.9 

-0.8 

0.0000 

2  0 

18397.1 

1636.4 

16760.7 

0.7162 

3  0 

17905.0 

2043.8 

15861.2 

0.6778 

4  0 

17049.3 

2765.5 

14283.8 

0.6104 

5  0 

13118.3 

6326.7 

6791.6 

0.2902 

Snow  depth=  0.0 

in  Density=  0.2 

g/cm''3 

G  4 

2340.1 

1228.6 

1111.5 

0.0475 

Snow  depth=  1 . 0 

in  Density=  0.2 

g/cm''3 

7  4 

2340.1 

1228.6 

1111.5 

0.0475 

Snow  depth=  10.0 

in  Density=  0.2 

g/cm^l 

8  4 

3763.6 

2314.2 

1449.4 

0.0619 

Snow  depth=  1 . 0 

in  Density=  0.2 

g/cm''3 

9  2 

9921.9 

1388.0 

8533.8 

0.3647 

Snow  depth=  10.0 

in  Density=  0.2 

g/cm''3 

10  2 

9921.9 

2314.2 

7607.7 

0.3251 

Snow  depth=  20.0 

in  Density=  0.2 

g/cm*3 

11  2 

9921.9 

3162.4 

6759.4 

0.2889 

Snow  depth=  10.0 

in  Density=  0.1 

g/cm*3 

12  2 

9921.9 

2003.0 

7918.9 

0.3384 

Snow  depth=  10.0 

in  Density=  0.3 

g/cm^3 

13  2 

9921.9 

2314.2 

7607.7 

0.3251 

Snow  depth=  10.0 

in  Density=  0.4 

g/cm^3 

14  2 

9921.9 

2003.0 

7918.9 

0.3384 

Snow  depth=  10.0 

in  Density=  0.2 

g/cm^3 

15  1 

9921.9 

1963.2 

7958.7 

0.3401 

16  3 

9921.9 

2407.2 

7514.6 

0.3211 

17  3 

9921.9 

2479.9 

7442.0 

0.3180 

18  3 

9921.9 

2910.1 

7011.8 

0.2996 

19  2 

9921.9 

2314.2 

7607.7 

0.3251 

20  2 

9921.9 

2314.2 

7607.7 

0.3251 

21  2 

9921.9 

2314.2 

7607.7 

0.3251 

22  2 

9921.9 

2314.2 

7607.7 

0.3251 

23  3 

9921.9 

2479.9 

7442.0 

0.3180 

24  3 

9921.9 

2910.1 

7011.8 

0.2996 

25  3 

9921.9 

2479.9 

7442.0 

0.3180 

26  3 

9921.9 

2479.9 

7442.0 

0.3180 

27  3 

9921.9 

2479.9 

7442.0 

0.3180 

28  3 

9921.9 

2479.9 

7442.0 

0.3180 

Snow  depth=  0 . 0 

in  Den3ity=  0.2 

g/cm^3 

29  0 

18397.1 

1636.4 

16760.7 

0.7162 

30  0 

18397.1 

1636.4 

16760.7 

0.7162 

31  0 

18397.1 

1636.4 

16760.7 

0.7162 

32  0 

18397.1 

1636.4 

16760.7 

0.7162 

33  0 

17049.3 

2765.5 

14283.8 

0.6104 

34  0 

8422.5  60120.7 

-51698.2 

-2.2092 

35  0 

17049.3 

2765.5 

14283.8 

0.6104 

36  0 

17049.3 

2765.5 

14283.8 

0.6104 

37  0 

17049.3 

2765.5 

14283.8 

0.6104 

38  0 

17049.3 

2765.5 

14283.8 

0.6104 

Tractn . 

Resist . 

[MPH] 

TFmull 

MRmull 

Slope= 

0 . 0  deg 

2340 

877 

29.8 

17950 

1636 

34.3 

17904 

2043 

20.9 

17049 

2765 

19.7 

13118 

6326 

8.2 

Slope= 

0 . 0  deg 

2340 

1228 

37.4 

Slope= 

0 . 0  deg 

2340 

1228 

37.4 

Slope= 

0 . 0  deg 

3763 

2314 

18.7 

Slope= 

0 . 0  deg 

9921 

1388 

39.0 

Slope= 

0 . 0  deg 

9921 

2314 

20.3 

Slope= 

0 . 0  deg 

9921 

3162 

16.9 

Slope= 

0 . 0  deg 

9921 

2002 

20.5 

Slope= 

0 . 0  deg 

9921 

2314 

20.3 

Slope = 

0 . 0  deg 

9921 

2002 

20.5 

Slope= 

0 . 0  deg 

2340 

1963 

16.5 

9921 

2407 

20.3 

9921 

2479 

20.2 

9921 

2910 

18.5 

9921 

2314 

20.3 

9921 

2314 

20.3 

9921 

2314 

20.3 

9921 

2314 

20.3 

9921 

2479 

20.2 

9921 

2910 

18.5 

9921 

2479 

20.2 

9921 

2479 

20.2 

9921 

2479 

20.2 

9921 

2479 

20.2 

Slope= 

0 . 0  deg 

17950 

1636 

34.3 

17950 

1636 

34.3 

17950 

1636 

34.3 

17950 

1636 

34.3 

17049 

2765 

19.7 

8422 

60120 

0.0 

0.642 

9.503 

17049 

2765 

19.7 

17049 

2765 

19.7 

17049 

2765 

19.7 

17049 

2765 

19.7 

66 


1994 


29-JUL-94  11:18:02  NRMMII  version  2.5.0  last  changed  12  May, 

TERRAIN  1  FILE=crrelcko.ter 

38  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  6  Ma 
SCENARIO  1  FILE=scenario.dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  2  FILE=vehicles/hmmwv.crl 

ID=  HMMWV,  M998,  4X4  FOR  CRREL  SNOW  MODEL  CHECKOUT  3  MAY  94 
Snow  model: CRREL 

Weight  [Lb] :  2®  1625  +  2®  2135  =  7520 

NTU/  - Gross-  (soil  only)  -  [Lb] - - Total-  [Lb]  -  Speed 

IBASE  Traction  Resist.  drawbar  DB-Coef.  Tractn.  Resist.  [MPH]  TFmull  MRmull  T 
Snow  depth=  10.0  in  Density=  0.1  g/cm''3  Slope=  0.0  deg 


1  0 

2924.6 

2924.5 

0.1 

0.0000 

752 

86 

29.8 

2  0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

3  0 

4462.7 

659.4 

3803.4 

0.5058 

4462 

659 

40.0 

4  0 

4237.4 

944 . 1 

3293.3 

0.4379 

4237 

944 

35.6 

5  0 

3142.0 

3189.6 

-47.6 

-0.0063 

3142 

3189 

0.0 

Snow  depth=  0 . 0 

in  Density=  0.2  g/cm^3 

Slope= 

0 . 0  deg 

6  4 

752.0 

131.6 

620.4 

0.0825 

752 

131 

40.0 

Snow  depth=  1 . 0 

in  Density=  0.2 

g/cm''3 

Slope= 

0 . 0  deg 

7  4 

752 . 0 

131.6 

620.4 

0.0825 

752 

131 

40 . 0 

Snow  depth=  10 . 0 

in  Density=  0.2 

g/cm^3 

Slope= 

0 . 0  deg 

8  4 

1243.7 

1245.5 

-1.8 

-0.0002 

1243 

1245 

0.0 

Snow  depth=  1 . 0 

in  Density=  0.2 

g/cm^3 

Slope= 

0 . 0  deg 

9  2 

2937.2 

510.4 

2426.8 

0.3227 

2937 

510 

40.0 

Snow  depth=  10.0 

in  Density=  0.2  g/cm^3 

Slope= 

0 . 0  deg 

10  2 

2937.2 

1245.5 

1691.7 

0.2250 

2937 

1245 

27 . 5 

Snow  depth=  20.0 

in  Density=  0.2 

g/cm''3 

Slope= 

0 . 0  deg 

11  2 

2937.2 

2504.4 

432.8 

0.0576 

2937 

2504 

10.3 

Snow  depth=  10.0 

in  Density=  0.1 

g/cm^3 

Slope= 

0 . 0  deg 

12  2 

2937.2 

813.1 

2124.1 

0.2825 

2937 

813 

40.0 

Snow  depth=  10.0 

in  Density=  0.3 

g/cm''3 

Slope= 

0 . 0  deg 

13  2 

2937.2 

1458.8 

1478.4 

0.1966 

2937 

1458 

26.4 

Snow  depth=  10.0 

in  Density=  0.4 

g/cm''3 

Slope= 

0 . 0  deg 

14  2 

2937.2 

1377.2 

1560.0 

0.2074 

2937 

1377 

26.9 

Snow  depth=  10.0 

in  Density=  0.2 

g/cm''3 

Slope= 

0 . 0  deg 

15  1 

2937.2 

1226.7 

1710.5 

0.2275 

752 

1226 

0.0 

16  3 

2937.2 

1255.2 

1682.0 

0.2237 

2937 

1255 

27.5 

17  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

18  3 

2937.2 

1499.2 

1437.9 

0.1912 

2937 

1499 

26.2 

19  2 

2937.2 

1245.5 

1691.7 

0.2250 

2937 

1245 

27.5 

20  2 

2937.2 

1245.5 

1691.7 

0.2250 

2937 

1245 

27.5 

21  2 

2937.2 

1245.5 

1691.7 

0.2250 

2937 

1245 

27.5 

22  2 

2937.2 

1245.5 

1691.7 

0.2250 

2937 

1245 

27.5 

23  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

24  3 

2937.2 

1499 . 2 

1437.9 

0.1912 

2937 

1499 

26.2 

25  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

26  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

27  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

28  3 

2937.2 

1304.5 

1632.7 

0.2171 

2937 

1304 

27.2 

Snow  depth=  0 . 0 

in  Density=  0.2 

g/cm^3 

Slope= 

0 . 0  deg 

29  0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

30  0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

31  0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

32  0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

33  0 

4237.4 

944.1 

3293.3 

0.4379 

4237 

944 

35.6 

34  0 

4237.4 

944.1 

3293.3 

0.4379 

4237 

944 

35.6 

35  0 

4237.4 

944.1 

3293.3 

0.4379 

4237 

944 

35.6 

36  0 

4237.4 

944 . 1 

3293.3 

0.4379 

4237 

944 

35.6 

37  0 

4237.4 

944 . 1 

3293.3 

0.4379 

4237 

944 

35.6 

38  0 

4237.4 

944 . 1 

3293.3 

0.4379 

4237 

944 

35.6 

67 


29-JUL-94  11:18:06  NRMMII  version  2.5.0  last  changed  12  May,  1994 
TERRAIN  1  FILE=crrelcko2 .ter 

54  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  6  Ma 
SCENARIO  1  FIIiE=scenario.dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  1  FILE=vehicles/mll3 . crl 

ID=  M113  FOR  CRREL  MODEL  CHECKOUT  3  MAY  94 


Snow  model : CRREL 


Weight 

[Lb]  : 

2®  11700 

=  23401 

NTU/ 

only) - [Lb] - 

- Total- [Lb] - 

Speed 

TRASE 

Traction 

Resist . 

drawbar 

DB-Coef . 

Tractn. 

Resist . 

[MPH] 

TFmull 

MRmull 

.Snow  deDth=  0 . 

0  in  Density=  0.2 

g/cm^3 

Slope= 

0 . 0  deg 

39 

0 

9360.1 

9360.9 

-0 . 8 

0.0000 

9360 

9360 

0.0 

40 

0 

18397.1 

1636.4 

16760.7 

0.7162 

17950 

1636 

34.3 

41 

0 

18397.1 

1636.4 

16760.7 

0.7162 

17950 

1636 

34.3 

42 

0 

-2521 . 7 

23082.5 

-25604.1 

-1.0941 

-2521 

23082 

0.0 

43 

0 

18397.1 

1636.4 

16760.7 

0.7162 

17950 

1636 

34.3 

44 

0 

18397.1 

1636.4 

16760.7 

0 . 7162 

17950 

1636 

34.3 

45 

0 

12125.3 

19540.4 

-7415.1 

-0.3169 

12125 

19540 

0 . 0 

0.711 

7.066 

46 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0 . 0 

0.711 

7.066 

47 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0.0 

0.711 

7.066 

48 

0 

8947 . 0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0.0 

0.711 

7.066 

49 

0 

12747.3 

11302 . 5 

1444.8 

0.0617 

12747 

11302 

1.7 

0.711 

7.066 

50 

0 

12913.7 

12956 . 0 

-42.3 

-0 . 0018 

12913 

12955 

0 . 0 

0.711 

7.066 

51 

0 

18158.0 

1833 . 6 

16324.4 

0 . 6976 

12045 

1833 

23 . 0 

52 

0 

12747 . 3 

11302.5 

1444.8 

0 . 0617 

12747 

11302 

1.7 

0.711 

7.066 

53 

0 

18158 . 0 

1833 . 6 

16324.4 

0.6976 

12045 

1833 

23.0 

54 

0 

18158 . 0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22.9 

55 

0 

18158 . 0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22 . 9 

56 

0 

18158 . 0 

1833.6 

16324.4 

0.6976 

12045 

1833 

23 . 0 

57 

0 

18158 . 0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22.9 

58 

0 

0.02316699.0********* 

-99.0000 

2340 

2316699 

0.0 

59 

0 

0.02316699.0********* 

-99.0000 

2340 

2316699 

0.0 

60 

0 

12747.3 

11302 . 5 

1444.8 

0.0617 

12747 

11302 

1.7 

0.711 

7.066 

61 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0 . 0 

0.711 

7.066 

62 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0 . 0 

0.711 

7.066 

63 

0 

12913.7 

12956.0 

-42.3 

-0.0018 

12913 

12955 

0 . 0 

0.711 

7.066 

64 

0 

12747.3 

11302.5 

1444 . 8 

0.0617 

12747 

11302 

1.7 

0.711 

7.066 

65 

0 

12296 . 6 

10346.2 

1950.4 

0.0833 

13817 

10346 

2.5 

0.711 

7.066 

66 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0.0 

0.711 

7.066 

67 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0 . 0 

0.711 

7.066 

68 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0 . 0 

0.711 

7.066 

69 

0 

12747.3 

11302.5 

1444 . 8 

0.0617 

12747 

11302 

1.7 

0.711 

7.066 

70 

0 

12913.7 

12956.0 

-42.3 

-0.0018 

12913 

12955 

0 . 0 

0.711 

7.066 

71 

0 

18158.0 

1833.6 

16324.4 

0.6976 

12045 

1833 

23 . 0 

72 

0 

12747.3 

11302.5 

1444.8 

0.0617 

12747 

11302 

1.7 

0.711 

7.066 

73 

0 

18158.0 

1833.6 

16324.4 

0.6976 

12045 

1833 

23 . 0 

74 

0 

18158.0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22.9 

75 

0 

18158.0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22.9 

76 

0 

18158.0 

1833.6 

16324.4 

0.6976 

12045 

1833 

23.0 

77 

0 

18158.0 

1833.6 

16324.4 

0.6976 

9217 

1833 

22.9 

78 

0 

0.02316699.0********* 

-99.0000 

2340 

2316699 

0.0 

79 

0 

0.02316699.0********* 

-99.0000 

2340 

2316699 

0.0 

80 

0 

12747.3 

11302.5 

1444 . 8 

0.0617 

12747 

11302 

1.7 

0 . 711 

7.066 

81 

0 

8947.0 

23975.3 

-15028 . 3 

-0.6422 

9910 

23975 

0 . 0 

0 . 711 

7.066 

82 

0 

8947.0 

23975.3 

-15028.3 

-0.6422 

9910 

23975 

0.0 

0 . 711 

7.066 

83 

0 

12913.7 

12956.0 

-42.3 

-0 . 0018 

12913 

12955 

0.0 

0.711 

7.066 

84 

0 

12747.3 

11302.5 

1444 . 8 

0.0617 

12747 

11302 

1 . 7 

0.711 

7.066 

85 

0 

17290.2 

1464.3 

15826 . 0 

0.6763 

17950 

1464 

38.2 

86 

0 

17049.3 

2765.5 

14283.8 

0.6104 

17049 

2765 

19.7 

87 

0 

17049 . 3 

2765.5 

14283 . 8 

0 . 6104 

17049 

2765 

19 . 7 

88 

0 

10946 . 5 

26279.6 

-15333 . 2 

-0.6552 

10946 

26279 

0 . 0 

0 . 642 

9 . 503 

89 

0 

10946 . 5 

26279 . 6 

-15333.2 

-0 . 6552 

10946 

26279 

0.0 

0.642 

9.503 

90 

0 

-2521 . 7 

23082 . 5 

-25604.1 

-1 . 0941 

-2521 

23082 

0.0 

91 

0 

13118.3 

6326.7 

6791.6 

0 . 2902 

13118 

6326 

8.2 

92 

0 

18158 . 0 

1833.6 

16324.4 

0.6976 

17950 

1833 

22.8 

68 


29-JUL-94  11:18:57  NRMMII  version  2.5.0  last  changed  12  May,  1994 
TERRAIN  1  FILE=crrelcko2 . ter 

54  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  6  Ma 
SCENARIO  1  FILE=scenario . dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  2  FILE=vehicles/hmmwv.crl 

ID=  HMMWV,  M998,  4X4  FOR  CRREL  SNOW  MODEL  CHECKOUT  3  MAY  94 
Snow  model : CRREL 

Weight  [Lb] :  2®  1625  +  2®  2135  =  7520 

NTU/  - Gross- (soil  only)  -  [Lb] - - Total- [Lb]  -  Speed 


IBASE 

Traction 

Resist . 

drawbar 

DB-Coef . 

Tractn. 

Resist . 

[MPH] 

TFmull 

MRmull 

Snow  depth=  0 , 

.  0  in  Density=  0.2 

g/cm^3 

Slope= 

0 . 0  deg 

39 

0 

2924 . 6 

2924.5 

0.1 

0.0000 

2924 

2924 

0.0 

40 

0 

2924.6 

2924.5 

0.1 

0.0000 

2924 

2924 

0.0 

41 

0 

4554 . 9 

481.7 

4073.3 

0.5417 

4554 

481 

40.0 

42 

0 

1616.7 

9306.3 

-7689.6 

-1.0226 

1616 

9306 

0.0 

43 

0 

3142.0 

3189.6 

-47.6 

-0.0063 

3142 

3189 

0.0 

44 

0 

4554.9 

481.7 

4073.3 

0.5417 

4554 

481 

40 . 0 

45 

0 

3013.6 

6670 . 6 

-3657.0 

-0.4863 

3013 

6670 

0.0 

0.711 

7.066 

46 

0 

1458.7 

2021.2 

-562.5 

-0.0748 

1458 

2021 

0.0 

0.711 

7.066 

47 

0 

1458.7 

2021.2 

-562.5 

-0.0748 

1458 

2021 

0.0 

0 . 711 

7.066 

48 

0 

1458.7 

2021.2 

-562.5 

-0.0748 

1458 

2021 

0 . 0 

0 . 711 

7.066 

49 

0 

3648.0 

3116.7 

531.3 

0.0707 

3647 

3116 
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29-JUL-94  11:18:12  NRMMII  version  2.5.0  last  changed  12  May,  1994 
TERRAIN  1  FILE=crrelcko3 .ter 

20  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  22,  J 
SCENARIO  1  FILE=scenario.dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  1  FILE=vehicles/mll3 .crl 

ID=  M113  FOR  CRREL  MODEL  CHECKOUT  3  MAY  94 
Snow  model : CRREL 

Weight  [Lb] :  2®  11700  =  23401 
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29-JUL-94  11:18:39  NRMMII  version  2.5.0  last  changed  12  May,  1994 
TERRAIN  1  FILE=crrelcko3.ter 

20  UNITS,  ID=Checkout  terrain  for  CRREL  winter  effects  model  22  J 
SCENARIO  1  FILE=scenario.dat 

"CRRELSNOW  "  ID=Dry, Snow, January  (CRREL  model) 

VEHICLE  2  FILE=vehicles/hmmwv.crl 

ID=  HMMWV,  M998,  4X4  FOR  CRREL  SNOW  MODEL  CHECKOUT  3  MAY  94 
Snow  model :  CRREL 

Weight  [Lb] :  2®  1625  +  2®  2135  =  7520 
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